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ABSTRACT 


This dissertation deals with an analysis of fatigue behavior at 
stress-concentration sites under cyclic loading. Specifically, the 
study treats sheet specimens containing a central hole under remote 
uniaxial cyclic loading. This investigation was conducted in 
three phases: an analysis of local cyclic stress-strain conditions 

at the .stress-concentration site, an investigation of fatigue crack 
initiation, and a study of the initial stage of crack propagation for 
cracks growing from the stress-concentration site. 

Because fatigue cracks in structures originate at stress 
concentrations, rational methods for estimating crack initiation 
must be based upon cyclic stress-strain conditions at typical stress- 
concentration sites. For the specimen configuration considered in 
the present study, this point of maximum stress and strain was on the 
boundary of the hole and consequently was in a uniaxial stress state. 
Cyclic stress and strain were estimated for this point by cycling 
uniaxial coupons between calculated limits for the local stress and 
strain. The adequacy of this "analytical-control" procedure was 
demonstrated by a comparison of calculated results with experimental 
results obtained on a companion specimen. Comparisons were made for 
constant-amplitude, two-level, and random loading. 


The first level of the two-level loading was designed to- produce 
residual stresses that altered local fatigue damage accumulation for 
the .second loading level. This nominal stress interaction was- typical 
of those experienced in structures. To obtain estimates of crack 
initiation that accounted for this stress interaction, the analytical- 
control tests, from the first phase were continued until, the uniaxial 
'coupons failed by fatigue. The accuracy of this procedure was 
illustrated by comparison of estimated and observed crack- initial ion 
periods. 

To complete this study of the role of stress concentrations in 
fatigue, the initial stage of ' crack propagation from a stress- 
concentration site was. investigated. Short cracks emanating from 
stress-concentration zones are influenced by the local stress -field 
corresponding to the stress concentration. Through the use of a 
Westergaard stress function for concentrated loads on a crack surface 
together with a boundary-collocation technique, a general procedure 
was developed for calculating stress- intensity factors for cracks 
growing from a hole in a sheet specimen. Stress -intensity factors 
were used together with a "master curve" characterizing crack growth 
for the sheet material to calculate crack-growth curves for cracks 
growing from the hole-. These calculations were made for local stress 
fields corresponding to the. constant-amplitude and two-level loading 
of the first two. phases of the study. The crack- initiation tests of 
the second phase were continued to obtain experimental crack 
propagation results. Eeasonable correlation was found between the 
calculated and experimental crack-growth curves. 


The generally close correlation "between calculated and observed 
results throughout this study demonstrated the utility of local stress 
analyses in the prediction of fatigue behavior for structural 
components containing stress concentrations. The procedures used 
herein are expected to be generally applicable for fatigue analyses 
involving more complicated configurations and loading than those 
considered in this investigation. 
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1. INTRODUCTION 


1. 1 General 

The importance of fatigue considerations in the design and use of 
structures has been well recognized for many years. However, in spite 
of intense efforts to develop analytical methods for predicting 
structural fatigue, most fatigue evaluations must be based upon ad hoc 
tests of prototype structures. These tests are usually time consuming 
and may result in expensive design changes. Analytical predictions of 
structural fatigue would be preferred at an early design stage of a 
new structure. 

One of the primary deterrents to the analysis of structural 
fatigue has been the problem of local stress analysis. Virtually all 
structures contain fasteners, fillets, or other such abrupt changes in 
configuration that act as stress raisers during service loading. The 
material at -each of these stress-concentration sites is highly stressed 
(often into the plastic range) during the cyclic loading and, con- 
sequently,- these sites become fatigue critical areas. Fatigue cracks 

originate at and propagate from these critical zones. To predict the 

} 

initiation of fatigue cracks, stress-strain conditions within these 

f 

zones must be determined. Moreover, during the early stages of crack 
propagation, a crack is under the influence of the local stress field. 
Consequently, the problem of local stress analysis may also complicate • 
the analysis of fatigue crack propagation. 
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1.2 Specific Objectives 

The present study was undertaken to investigate typical effects of 
stress concentrations on the fatigue behavior of sheet specimens. This 
study was conducted in three phases: a study of local cyclic stress- 

strain conditions at the stress-concentration site, an investigation 
of crack initiation, and an analysis of the initial stage of crack 
propagation at this site. 

Because fatigue cracks originate at the most highly stressed zone 
of a stress-concentration site, attention was restricted to this zone 
during the first phase of the study. As a 'further simplification, thin 
specimens containing a central hole with inplane loading were con- . 
sidered. For this type of specimen, the vicinity of the stress 
raiser is in a state of plane stress and the fatigue critical zone on 
the stress-free boundary of the hole is under simple uniaxial stress. 
Conditions at this critical point were referred to as "local" 
conditions. As a result of this uniaxial stress state, local 
conditions were investigated by simulating local cyclic stress- strain 
behavior in unnotched uniaxial specimens. This simulation was 
performed by cycling the unnotched specimens between prescribed 
local elastoplastic stress- strain limits. An equation relating local 
stress- and strain-concentration factors for nonlinear material 
behavior was used to calculate these limits. For an evaluation of 
this "analytical-control" procedure, local cyclic stress-strain 
conditions were also determined experimentally using a "companion- 
specimen" method. The calculated and experimental results were compared 
for constant and variable -amplitude fatigue loading. 
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To demonstrate the importance of analyzing fatigue crack initia- 
tion from a local cyclic stress-strain approach,, loading levels and 
sequences were selected to produce stress interactions. These stress 
interactions were caused by residual stresses and damage acceleration 
at the stress-concentration site and were typical of those found in 
service. Fatigue crack- initiation periods for the notched specimens 
were estimated by continuing the -analytical -control tests fr.om the 
first phase until cracks appeared in the unnotched specimens. The 
adequacy of this prediction method was assessed by comparing estimated 
results and observed crack-initiation periods in notched specimens. 

The effect of the stress concentration on fatigue exists beyond 
the initiation of a fatigue crack at the stress-concentration site. 

r * 1 

The high stresses in the immediate vicinity of a stress raiser 
increase the crack propagation rates for small cracks originating at 
the critical site. In addition, for variable -amplitude loading 

■i * ' ^ 

involving local plasticity, residual stresses may exist in the 
vicinity of the stress concentration that alter crack-growth behavior. 
Consequently, a third phase .-of this study investigated the effects of 
stress raisers during the initial' stage of crack propagation. 

i 

Through the use of a Westergaard stress function for concentrated 

J < 

loads on a free crack together with a boundary- collocation procedure, 
a general procedure was developed for calculating stress- intensity 
factors for cracks growing from a central hole in sheet specimens. 
Stress- intensity factors for cracks propagating through residual 
stress fields were determined with the aid of elastic and approximate 
elastoplastic analyses of stresses near the circular hole. These stress 
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intensity factors were used together with a "master curve," character- 
izing the crack- growth behavior for the specimen material, to calculate 
cr.ack-growth curves for several - types of constant and variable- 
amplitude loading considered in the first two phases of this study. 

The accuracy of the estimates of stress-interaction effects on initial 
crack propagation behavior was* determined by comparing estimated 
crack-growth curves with test results. 

The correlation between calculated and experimental results 
demonstrated in this study illustrates the utility of local stress 
analyses for the prediction of fatigue lives for structural components 
containing stress concentrations. The procedures used in this study 
to estimate fatigue behavior are expected to be generally applicable 
in fatigue analyses involving more complicated configurations and 
loading sequences than those considered in this investigation. 


2. LITERATURE REVIEW 


2 .1 General 

Fatigue failures in machine and. structural components 
characteristically originate at stress concentrations such as bolt 
holes, fillets, or other abrupt changes in configuration. Because of 
this localized nature of fatigue crack initiation, rational methods 
for predicting fatigue must be based upon analyses of local stress 
and strain at the stress-concentration sites. In weight-critical 
structures, because of the high design stresses, local conditions 
often extend into the plastic range. This literature review will 
briefly survey the problem area of elastoplastic analysis of cyclic 
stress-strain conditions at stress cpncentrations. In addition, 
methods for predicting fatigue crack-initiation periods based on local 
stress are reviewed in the second part of this literature' survey. 

To complete the review of the effects of stress concentrations on 
fatigue, the third part of this review deals with the initial stage of 
crack propagation; for 1 cracks growing from stress-concentration sites. 

2 . 2 Local Cyclic Stress-Strain Studies 

The analysis of local elastoplastic stress- strain behavior at a 
stress-concentration site is quite complicated even for monotonic 
loading. For cyclic loading, the local stress and strain depend also 
upon the loading sequence which severely complicates theoretical 
analyses. Furthermore, many materials exhibit stress-strain properties 
that are altered by cycling. As a result of these complicating factors. 
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numerical or 'approximate techniques have been used to analyze local 
conditions for cyclic' loading. 

2.2.1 Numerical Analyse s 

* S 

1 ^ I 

Several numerical analyses have been conducted 'for plane specimens 
under cyclic loading. Denke used a finite-element approach to 
obtain a numerical solution for the elastoplastic stress distribution 
near a central hole in a plate during one cycle of loading. Material 
behavior was assumed to be elastoplastic during loading and entirely 
elastic during unloading. From this solution residual stress 
distributions were determined and used to estimate fatigue behavior. 

Ibrahim, McC allion, and Dudley (jf] also investigated the elasto- 
plastic deformation around a circular hole under one cycle of loading. 

A finite-difference approach together with a relaxation method were 
used to obtain a plane-strain solution. Reasonable correlation was 
shown between this solution and test results on a sheet containing a 
central hole in plane stress. The elastoplastic analysis was con- 
ducted for both loading and unloading to account for the reversed 
yielding associated with unloading. 

Isakson, Armen, and Pifko &] applied a finite-element analysis 
throughout three cycles’ of reversed loading of an edge-notched sheet 
specimen. Specimen configuration, material behavior, and loading 
levels were identical to those used by Crews [V] in an experimental 
study of local behavior and a close correlation was found between 
numerical and experimental results. 

A finite-element technique developed by Jordan [V] was applied by 
Mowbray and Slot \§ \ to investigate stress and strain redistribution 


t 
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in notched specimens during , cyclic loading. Mowbray and Slot selected 
the specimen dimensions, material, and loading used by Blatherwick and 
Olson [Y] in an experimental analysis' of strain- distribution. ^ Cyclic 
stress-strain diagrams were used as effective 'stress-strain curves to 
determine stress and strain redistributions resulting from the cyclic 
strain softening exhibited by the- specimen material. The numerical 
results were found to correlate well with measured strains reported by 
Blatherwick and Olson. 

2 .-2. 2 Stress- and Strain-Concentration, Studies 

Prom a fatigue- point of view, attention can be restricted to the 
most highly stressed zone of a stress-concentration site without 
seriously affecting the usefulness -of stress-concentration studies. 
This zone becomes the fatigue critical site of the specimen for cyclic 
loading. If fatigue 1 cracks develop, they will originate in this zone. 
In this study attention was consequently limited to this most critical 
location. Conditions at this -site were referred to as "local" 
conditions. Furthermore, for the case of plane stress considered in 
this investigation, the critical zone reduces to a point and the 
behavior at this location can be conveniently described in terms of 
stress- and strain-concentration factors. 

In the second part of this literature survey several fatigue 
prediction methods based on the use of stress- or strain-concentration 
factors are reviewed. All of these methods are based upon either the 
modified Stowell [Y] or the Weuber [9] equations for elastoplastic 
stress- and strain-concentration factors. 

1 r 

1 * ' 

» 

' l 

r 

4 V * * 

t . 
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Stowell investigated the problem of elastoplastic behavior of a 
sheet with a circular hole and with uniform loading at infinity. To 
account for local plasticity., Stowell introduced the ratio Eg/E^ 
into the elastic solution where E s is the secant modulus at the 
critical location and Em is the secant modulus corresponding to the 
uniform stress remote from the hole. This formulation led to the 
stress-concentration factor for elastoplastic behavior of 


K-. = 1 + 2 — 

U TTT 

lb CO 


( 1 ) 


The corresponding strain-concentration factor was found to be 


' Eg, 

K e = 1 + 2 
e E„ 


Em 

Es 


(2) 


Stowell compared his results with the experimental • data obtained by 
Griffith jjLoJ for wide 2024-T3 aluminum alloy sheets also with a 
circular hole under remote uniform tension. Excellent correlation 
was found for stress-concentration factors. However, theoretical 
strain-concentration factors were noticeably smaller than the experi- 
mental results. 

Hardrath and Ohman generalized the stress-concentration 

factor equation (l) -by presenting the following expression in terms 

of the' elastic stress-concentration factor Kj 

; < 

\ - 4 

'■ ■ , *' K i -f - 1 + - nf 

' <, ' ’ 

„ ' , V 

l * 

* . • ' ’ * ■ . I 


( 3 ) 
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where a . is the maximum local stress and S is the nominal stress 

, » 

at a remote location. This equation was used to calculate stress- 
concentration factors for sheet specimens containing notches and 
'fillets, and results were compared with experimental data for 2024-T3 
aluminum alloy. In addition, Hardrath and Ohman compared calculated 
results with experimental stress-concentration factors determined by 
Box [12], In general, equation (3) correlated very closely with 
experimental results. 

Crews [ 13 ] and Crews and Hardrath [iCJ considered the application 
of the generalized Stowell equation (3) for cyclic loading conditions. 
The nominal stress sequence was divided into monotonic excursions and 
equation ( 3 ) was applied for each excursion. The initial stress-strain 
state for each monotonic excursion was taken as the final state from 
the previous excursion. A generalization of the Stowell equation used 
in this procedure can be written as 



1 + (% - 


W 


where the subscript "i" refers to the ith monotonic nominal stress 
excursion ASq. Results from equation (it-) agreed closely with 
experimental results in [l33 und [l 1 ^ for 2024-T3 aluminum alloy 
specimens under one cycle of loading. 

Equation (it-) was rewritten by this author as 

agj aSj 

- (Bp - l)^i 


ACTi = 


(5) 
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where Aa^ and Ae^ are the independent variables. For given values 
of Kip and E this equation was used to establish the extreme values 
of local stress and strain for each nominal stress excursion A3j_. 

These extreme conditions for local behavior were used to control the 
cyclic loading of unnotched specimens. The local stress-strain curves 
found by this procedure were compared with experimental results for 
sheet aluminum and steel specimens with stress-concentration factors 
of 2, and 6 throughout one cycle of completely reversed loading. 

Neuber presents a general theory relating stress and strain- 
concentration factors for shear-strained prismatical bodies and 
arbitrary stress- strain laws. Starting from the observation that 

K e >K T .>K ff (6) 

Neuber postulated the existence of a function, which was neither stress 
or strain, but a certain combination of both with the characteristic 
that its concentration factor had the same value for all stress-strain 
laws and therefore was equal to K^. This function was called the 
, "leading function" E(t) and was defined by 


ff(r) 

N(t w ) 


= % 


( 7 ) 


where t and Tjj are local and nominal shear stresses, respectively. 

The leading function was determined from the assumption that the 
curvature of the notch root remained unchanged for both elastic and 
general stress-strain behavior. This approach resulted in 


where F(t) was a strain function. From equation (7) and. the 
definitions 


, r T max 

„ F ( T max) 

h e _ — 

e F( 

Neuber rewrote (8) as 

% = ’ ( 9 ) 

For an evaluation of the leading -function method leuber &€ applied 
equation, (8) to determine stress-concentration factors for sharp 
notches in pure shear and compared these results with an exact 
solution. Satisfactory agreement was shorn. 

Neuber extended the leading-function method to cases of plane 

i 

stress in [17] through the use of the deformation theory of Nadai jl8] 
to express the strain function for 1 shear by the strain function for 
tension. To demonstrate the accuracy of the leading function procedure 
for plane stress, Neuber applied this method to determine stress- 
concentration factors at a circular hole in a sheet subjected to. uniform 
radial tension. For comparison, Neuber followed the approach of 
Budiansky and Mangasarian [19] to obtain an exact solution for this 



problem corresponding to several levels of strain hardening. For all 
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cases, the values obtained by the leading-function method were very 
near to those obtained for the exact theory, demonstrating the accuracy 
of the leading-function procedure for plane - stress. 

The Weuber equation (9) was first applied in cyclic loading by 
Peterson |j2o] and later by Dolan 1 Jj2l)'. These authors as well as Manson 
and Hirschberg 0*3. Topper, Wetzel, and Morrow [ 25 ] and Wetzel (2^+j 
have used the Webber equation in fatigue studies of notched specimens. 
These papers will be discussed in the literature review of fatigue 
prediction methods, and are listed here to demonstrate the wide use 
of equation (9) in applications with cyclic, loading. 

Crews [ 15 ] rewrote equation ( 9 ) in terms of local stress and 
strain as 




(% ASj) 2 
Aei E 


( 10 } 


where Aa^ and Aeq are the local stress and strain excursions 
corresponding to the ith nominal stress excursion Equation (10) 

was used to calculate limiting conditions for each cycle of local, 
stress and strain, and unnotched specimens, were cycled between these 
limits to determine the local stress-strain curves. This procedure 
will be described in detail in Chapter 3» 

Dixon and Stannigan [ 25 ] and Dixon \z 6\ experimentally investi- 
gated, the elastoplastic conditions around crack tips for loading and 
( » 

unloading using a photoelastic coating technique. As a result of these 
studies, they constructed expressions for stress and strain around the 


■ 1 
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crack tip during the first load cycles. For the case of notched 
specimens 


Act = 



( 11 ) 


would be given, where Aa 0 is the stress excursion corresponding to 
completely elastic behavior. Since 


and 


Ao 0 = Kj A3 


E s - 


Act 

Ae 


equation (ll) can be rewritten as 

(Em AS ) 2 

Aa = 

Ae E 

Although developed by an entirely different approach, this expression 
is identical to equation (10), developed from the Neuber equation. 

The close correlation found in [ 25 ] and. [26] for results obtained by 
equation (ll) and by photoelastic methods further demonstrated the 
accuracy of the Ileuber equation. 

2 . 3 Methods for Predicting Crack Initiation in Notched Specimens 
2.3.I General 

Fatigue prediction methods for notched specimens can be classified 
into two general categories. The first of these is based upon constant 
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stress amplitude tests to generate a series of S-N curves for 
specimens with various stress-concentration factors. To estimate 
the fatigue "behavior for a given specimen configuration, the stress- 
concentration factor must he calculated or determined experimentally 
and a constant stress amplitude fatigue life is obtained from the 
family of S-N curves hy interpolating for the desired Kqi value. 

The second hasic approach is "based on the assumption that a stress- 
concentration site and an unnotched tensile specimen experience 
fatigue cracks in the name number of cycles, if they are each cycled 
through the same stress-strain history. Thus, once local cyclic 
stress-strain conditions are determined for a stress-concentration 
site, unnotched specimens can be used to estimate the fatigue crack- 
initiation period for the stress raiser. The fatigue behavior of 
unnotched specimens can be characterized by a family of S-N curves. 

The fatigue tests for the various Kjj levels required by the 
first approach may be quite expensive and time consuming in contrast to 
simple laboratory tests of unnotehed specimens used in the second 
approach. Furthermore, loading sequences in variable-amplitude tests 
affect fatigue behavior of notched specimens and, as a result, a 
family of S-N curves for various Kq, values, in general, apply only 
for a single loading sequence. On the other hand, the effects .of 
loading sequence could he incorporated in the analysis of local stress 
and a single family of S-N curves for unnotched specimens could he used 
for all types of loading. For these reasons the local stress analysis 
approach will be followed in this study and will he reviewed in the 
following sections. 
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Because of simplifying approximations that can he made for local 
conditions corresponding to tension-tension and completely reversed 
loading, the following fatigue prediction methods are grouped according 
to type of loading. 

2. J. 2 Prediction Methods for Tension-Tension Loading 

Gunn [ 27 J applied the modified Stowell equation (3) together with 

.a mon'otonic stress-strain curve to estimate the maximum stress at a 

stress raiser. Gunn assumed that the local stress range was elastic 

after the first cycle, that isj, Acr = AS. From the estimates of 

maximum local stress and Aa, Gunn predicted the manner in which 

yielding lowered the' local mean stress for cases of tens ion-tens ion 
> 

loading. By these estimates of local stress' behavior, Gunn predicted 

« i ** ' f 

alternating- mean fatigue diagrams for notched specimens from similar 
diagrams for unnotched specimens. Although no provisions were made to 
account for reversed yielding or variations in cyclic material behavior, 
qualitative agreement was found between predicted and experimental 
fatigue lives. 

Independently, Smith Jj 28 ] developed a similar method for consider- 
ing the effect of plasticity on life predictions. Smith assumed that 
the strain at the concentration site was proportional to nominal 
stress and calculated the strain from e = KpS/E. The corresponding 
notch-root stresses were found from a monotonic stress-strain curve. 
Residual stress was taken as the difference between the maximum local 
stress and the elastic local stress range, as in &?]• However, to 
produce agreement with life predictions based on these results, the 
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residual stresses were multiplied by an assumed ratio between 

i ■* 

plastically and elastically deformed material. Local stresses found 

by this procedure , were used together with an S-N curve for unnotched 
‘ ' V 

specimens to predict fatigue lives of notched specimens under several 
, • . ' 1 . ’ : ‘ * , • 

levels of reversed (R'= . 0 ) loading. Miner’s linear cumulative damage 
theory |j29} was used to combine these estimated lives for predictions 
of fatigue behavior under two-level loading. In addition. Smith [j 5 o] 
applied the linear strain theory from [28] to predict fatigue lives 
for notched specimens under repeated (R = 0 ) block- loading and 
demonstrated that the highest load level in the spectrum determines 
the magnitude of the residual stress to be used in estimating the 
fatigue lives at all other load levels. 

Ripp, Gomza, and Adee [31] presented a procedure similar to those 
of Gunn and Smith, however, the Heuber stress-concentration factor 
% 0*0 was use ^ instead of Krp. 


% 


1 + 


% " 1 



(12) 


where A' is a material constant, p is the notch-root radius, and w 
is the flank angle of the notch. Because K™- < Kj, the notch-root 
strains estimated by Ripp et al. were lower than those obtained by 
either Smith or Gunn. The comparisons of predicted and observed 
fatigue lives in [52] reflected the large errors in strain based upon 
equation ( 12 ), especially for the low-life range. 
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Each of the preceding methods is based upon the assumption of 
elastic cyclic behavior at the notch root for cycles beyond the first. 
For R > 0 and for low nominal stresses this assumption is justified. 
However, for R = 0 and moderately high nominal stresses, reversed 
yielding occurs upon unloading (see [ 15 J ) • Consequently,, the fatigue 
prediction procedures for tension-tension loading must be restricted 
to low nominal stress ranges for best results. 

2.5-3 Prediction Methods for Completely Reversed Loading 

Peterson [2(0 modified the generalized Stowell equation ( 3 ) for 
the fatigue problem by replacing Kj by Kj 

K ff = 1 + (Kf - 1)^ (13) 

where Kp is the fatigue notch factor, obtained as the ratio of 
fatigue limits for unnotched and notched specimens. Peterson used 
equation ( 13 ) together with a cyclic stress-strain curve to estimate 
local stresses corresponding to specific levels of reversed (R = -1) 
nominal loading. Starting' with an S-N curve for unnotched specimens 
■under reversed loading, Peterson constructed an S-N curve for notched 
specimens in the following manner. The fatigue life was found for 
each calculated a (where cr is local stress) from the S-N curve 
for unnotched specimens and the S used to calculate 0 was plotted 
at this new fatigue life. The procedure was repeated for a range of 
S to obtain the S-N curve for notched specimens. The procedure was 
applied for edge-notched sheet specimens of SAE TlJO (normalized) 
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steel with elastic stress-concentration factors of 2 and 1)-. Close 
* * 

correlation was found "between calculated diagrams and the observed 
S-N curves in - Qjj] • 

Dolan [2X] suggested a procedure for estimating fatigue lives for 
notched members based on the Neuber equation ( 9 )* Dolan rewrote 
equation (9) in terms of cyclic amplitude of elastic strain as 


Krp AS 
E 




(I*) 


where Ae 0 and Ae^ are, respectively, the elastic and total ranges 
of local -strain and F - \J Ae e /Ae-f For given Ej and AS the left 
side of equation (l^) can be evaluated. The right side of equa- 
tion (ik) can be found by successive approximations from a hysteresis 
loop for the material. The estimate of fatigue life corresponding to 
Ae^ can ‘be found from a strain-life fatigue diagram for unnotched 
specimens. Dolan did not investigate the accuracy of this procedure 



Manson and Hirschberg (22]] presented an approximate analysis to 
estimate the number of cycles required to initiate a crack and the 
number of cycles required to propagate the crack, to failure for 
notched specimens. Cracks were assumed to occur at the root of a 
notch after a certain number of cycles dependent only on the localized 
surface strain as in the previously described methods. After the 
crack reached an "engineering size" it was assumed to propagate in a 
manner similar to those in unnotched specimens. The number of cycles 
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to initiate a crack was designated Ef 0 and was calculated for unnotched 
specimens using 

(15) 

from [31 ] , -where Wf is the total life to fracture. Manson and 
Hirschberg estimated the local strain range using the Stowell equation 
in the form suggested by Peterson equation (13) and Neuber's equa- 
tion ( 9 ) after replacing by K^. 

Kf = \/KA ( 16 ) 

Local strain estimates were used for establishing predictions of crack 
initiation at stress concentrations. Life estimates for crack 
propagation were also made with the assumption that crack growth 
depended only on the nominal strain range and not on the strain- 
concentration factor for the notch. Total fatigue life of a notched 
specimen was the sum of the- crack-initiation period and the propagation 
period. To evaluate this procedure Manson and Hirschberg calculated 
crack- initiation and propagation periods using equation (l 6) and compared 
their results with the observed behavior of 7075 -T 6 aluminum alloy and 
AISI 7130 steel specimens having stress-concentration factors of 2 and 3* 
Close correlation was found between calculated and observed results. 

Topper, Wetzel, and Morrow Q>3j wrote equation (l 6 ) in the form 


K^AS As E) 1 / 2 = (Acr Ae E) 1 ^ 2 


(17) 


( 4 T 


i 


where As is the nominal strain range. Equation (17) relates the 
nominal cyclic stress-strain behavior of notched specimens to the 
local cyclic stress-strain behavior. For completely reversed constant- 
amplitude tests detectable cracks "were assumed to form in notched and 
unnotched specimens in the same number of cycles provided 

l /2 1/2 

Kp (AS Ae E)- for the notched specimens was equal to (A cs Ae E') 

for the unnotched specimens. To apply this assumption. Topper et al. , 
used unnotched specimens to produce a "master curve" relating 
(Act As E) 1 ^ and fatigue' lives for unnotched specimens. The life 
(initiation phase) can be predicted for a notched specimen by entering , 
the master plot at the value of (Act Ae E)"^ 2 equal to Krji(AS Ae E) 1 ^ 2 
for the notched specimen test condition. Life predictions were made by • 
this procedure and compared with observed fatigue lives for edge-notched 
7075-T6 aluminum alloy specimens reported in QjjJ. The close correlation 
indicates that for this material the crack propagation phase can be 
neglected. 

All four of the preceding methods apply, in general, only to 

constant-amplitude completely reversed loading. For this special type 

of loading, local stresses and’ strains can also be assumed (see jjL53) 

to experience complete reversals and, as a result, local mean stresses 

are zero. Because of this simplification, maximum and minimum values 
* 

of local stress and strain can be calculated directly from a cyclic 
stress-strain curve (generated also by reversed loading) by either the 
modified j Stowell or the Heuber equation. Although not mentioned in 
this assumption was’ implicit in the procedure suggested by Dolan. 
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2.3.* Prediction Methods for Arbitrary Loading 

f * 

To predict fatigue of notched specimens under arbitrary types of 
loading using the local stress analysis approach, local cyclic stresses 
and strains must be determined, in general, for each loading cycle. 

This procedure is complicated by the inherent history dependence of local 
plasticity and by cyclic variations in stress-strain properties. 

Efforts to simplify this complex problem have led to the approximate 
analyses and experimental .studies reviewed in this section. 

The Ranalli method was described in [35] by Silverman, Hooson, 
and Saleme. This procedure is an extension of Peterson's method [20] 
to the case of nonfully reversed loading and is based upon semi- empirical 
methods for determining the strain amplitude and mean stress at the notch 
root. An equivalent fully reversed strain amplitude giving the same 
fatigue damage was determined and a fatigue life estimate was obtained 
for the equivalent strain amplitude from fully reversed constant strain 
amplitude fatigue curves. 

The Rana ll i method was applied to estimate fatigue lives for 
constant -amplitude tests and for simple spectrum loading. The linear 
damage theory (Miner's Rule) was used for the spectrum loading tests 
to account for fatigue damage at each stress level. Although based upon 
semi-empirical methods for analyzing local behavior, the' Ranalli method 
was- generally concluded to he of considerable value as a practical 
fatigue design method. 

Wetzel [2^t| applied equation (9) to calculate control conditions 
under which unnotched specimens were made to simulate fatigue behavior 
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of material at a notch foot. This procedure was similar to that used 

1 i 

by Crews _ { 15 } . However, Wetzel replaced Ej by Kf in the Weuber 

4 * ' 

equation to calculate the control conditions. This substitution improves 
1 the fatigue-prediction results for const ant -amplitude loading, but 

1 1 i 

* * 

could result in large errors for local stress-strain calculations, 
especially for high • values where K^, and Kg. are usually quite 

j- f 

- 1 * / « ' 

different.- -Endpoint conditions for each half-cycle of local stress- 
strain behavior were calculated from equation ( 9 ) and an unnotched 
specimen was suscessively cycled between the limits. This procedure 
can be applied for complicated types of loading, since each half-cycle 
of loading is handled separately. Wetzel applied the procedure to 
simulate local fatigue behavior for R = 0 loading, and close correla- 
tion was found between calculated and observed results. 

In [l], local stress conditions were determined by a numerical 
elastoplastic analysis but, in addition, experimental methods were used 
to find local stresses for const ant -amplitude R = 0 loading. 

Unnotched specimens were cycled between these experimental' limits and 
the resulting fatigue lives, were compared with lives found from fatigue 
tests of notched specimens. For the two test conditions considered, 
close correlation was found. 

A more extensive experimental study of local stress behavior was 
presented in [«]• Local strains were measured while notched specimens 
were loaded cyclically arid these measured strains were reproduced in 
tests of unnotched specimens to detemnine the corresponding local stress 
histories. For the 2021— T3 aluminum alloy investigated, local stresses 
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were found to reach stabilized levels. Fatigue lives for notched 
specimens were estimated from these stabilized local stresses and 
fatigue data for unnotched specimens. 

2. 4- Analysis of Crack Propagation Hear Stress Concentrations 

Crack initiation may occur in the early stages of fatigue life with 
the result that the crack propagation phase must be considered for 
accurate estimates of total fatigue life. Fatigue cracks initiate at 
areas of high stress concentration and grow away from these areas as 
cyclic loading is continued. The initial stage of this crack propagation 
is influenced by the stress field due to the stress raiser, but as the 
crack grows it gradu ally escapes this influence and propagates solely 
under the influence of the applied nominal stress. The behavior of 
fatigue cracks under . applied nominal stress has received much attention 
in the literature (for example, jj)6], and [ 38 ]) and, as a result, 

was not dealt with here. On the other hand, fatigue crack growth near 
stress concentrations has received' only limited attention. 

. Manson and Hirschberg (22^ discussed the transition period between 
crack initiation .and crack propagation under nominal stress. This 
period was included in their analysis by defining N 0 , not as the number 
of cycles to form a microcrack, but rather as the number of cycles to 

i 

form an engineering size crack according to the equation (15) • While 
Manson and Hirschberg recognized that this simplification was highly 
questionably, the results associated with this' assumption were compared 
with observed behavior to assess its validity. In general, the 

t ■ 

comparisons made in [ 22 ]] supported their simplifying assumption, 
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however, the validity of this approach is contingent upon the appropriate 
selection of the engineering size crack. 

Because this transition phase of crack growth depends upon the 
stress distribution around the stress concentration, loading sequences 
that alter this stress distribution through the creation of residual 
stresses also influence this initial phase of crack growth. While 
several .studies have been conducted to establish the effects of prior 
loading on fatigue crack initiation or total fatigue life (for example, 
[39] and Qs- 9 ^ ) , the effects of this loading on crack propagation have 
,not been analyzed. Taira and Murakami M measured crack growth from 
V notches in plates that had been statically loaded in tension prior 
to fatigue testing. 

. Prior tensile loading was found to suppress the initial stage of 
crack growth. This decrease in crack growth rate was attributed to 
compressive residual stresses produced by the prior loading. However, 
no analysis of this effect was presented. 

In summary, although several investigators have mentioned the 
transition phase of crack growth from a stress raiser, no deliberate 
study has been made of this phase of fatigue, for some applications 
this phase may represent a significant part of the fatigue life, and, 
as a result, should be considered in fatigue predictions. 



3. ANALYSIS OF LOCAL CYCLIC STRESS AND STRAIN 


3.1 General 

The present study of local stress and strain was undertaken to 
demonstrate a procedure for estimating local conditions and to 
establish the adequacy of this procedure by comparison of calculated 
and experimental results. In addition, the local stress- strain 
results developed in this study were used in the next chapter, -which 
deals with crack- initiation behavior at stress-concentration sites. 

The present interest in the analysis of local stress and strain 
stems from the fact that fatigue is a very localized process. Without 
adequate procedures for calculating the local conditions at the 
fatigue critical sites of a structure, predictions of fatigue life 
are impossible. For weight critical structures, the analysis of local 
cyclic behavior is quite complicated. The history dependence of local 
plasticity combined with the cyclic strain hardening or’ softening 
experienced by most materials virtually preclude exact solutions for 
local stress and strain beyond the first few loading cycles. 

For the purposes of this study the problem of local stress 
analysis was simplified by considering the specimen configuration shown 
in figure 3»1* (This specimen will be referred to herein as the 
"notched" specimen. ) Because this specimen was thin and was stressed 
by inplahe loading,' a plane stress condition’ existed. The elastic 
stress distribution [I 2 ] oh the transverse axis is shown in figure 3.2 
.for ^the specimen dimensions given in figure 3* 3* The fatigue critical 

, 4 1 
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points for this ■ configuration lie on the boundary of the hole at 
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Normalized ' stress 
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9=0 and a. Because these points were on the stress-free boundary 
of the hole, they were in a uniaxial state of stress which leads to 
significant simplifications in this study of cyclic behavior. 

Throughout the remainder of this chapter, attention will be restricted 
to these points and conditions at these critical points will* be 
referred to as "local" conditions. Based on net-section nominal stress, 
the elastic stress-concentration factor for these points is 2.57 M- 
In addition, the comparison in figure 5*2 illustrates no effect on 
local conditions by the finite .specimen width shown in figure 3*2. 

As a result of the uniaxial stress state at the local points, the 
. * 

local cyclic conditions can be simulated in simple unnotched uniaxial 
•specimens. ■ Local behavior was determined using this procedure. 

Two methods were used to control the simulation. In the first method, 
referred, to here as analytical-control method, limiting conditions 
involving local stress and strain were calculated and the unnotched 
uniaxial specimens were cycled between these limits. The second method 
: was entirely experimental in nature, and was used to establish- the 
accuracy of the analytical- control method. This experimental method, 
called the companion- specimen method, was based upon the experimental 
determination of local cyclic strain limits and the cyclic testing of 
unnotched companion specimens between these limits. 

Two types of cyclic loading were considered in this investigation. 
To demonstrate the variations in local stress and strain due to cyclic 
strain hardening and to establish the ability of the analytical- control 
method to account for this material behavior, tests and analysis were 
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conducted for constant-amplitude cyclic loading. The effects of 
loading sequences on local behavior were studied for several simple 
cases of variable-amplitude loading. 

i 

Details of procedures and apparatus are presented in the following 
sections. Discussions and comparisons are made between the calculated 
and experimental results for an evaluation of the analytical-control 
procedure. 

3.2 Analytical-Control Method 

The procedure described in this section is based upon the simula- 
tion of local behavior in unnotched specimens. Because the fatigue 
critical points for the notched specimen configuration considered in 
this study experienced only uniaxial stress, unnotched specimens under 
simple uniaxial loading were used for the simulation of local stress 
and strain. These plain specimens were cycled between calculated 
limits , involving local stress and strain, to perform the simulation. 

Two equations are currently available for calculating the limiting 
conditions for local stress and strain. These are the generalized 
Stowell and the TTeuber equations. 

The Stowell relationship given in the literature review as 
equation (4) 


i±S3A 



was rewritten in terms of Aoq and Ae^ as 


( 4 ) 
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Ae. 7 ASj 




( 5 ) 


Similarly the Neuber equation was written as 


> * * 




(K T A5j)g ’ 
E Ae^ 


( 10 ) 
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For the ^purpose ,of establishing the extreme values of local stress and 
strain for the analytical- control method, -either equation (5) or 
equation (10) -could have been used. However, based upon the simplicity 
of equation (10) and upon the evaluation of the analytical- control 
procedure for constant-amplitude loading performed by the author in 
[15], equation (10) was selected for this study. 

The cycle-by-.cycle procedure for applying the analytical-control 
method can be explained with the aid of figure 3-b. An arbitrary 
nominal stress cycle has been divided into monotonic excursions in 
figure 3-b(&), and the corresponding estimates for the local stress- 
strain behavior are presented in figure 3<Mb). To illustrate the 
procedure, equation (10 ) has been evaluated for a given Kj and E 
and plotted in the first quadrant of figure 3-b(jo) for ASj_. This 
curve for equation (10) represents the locus of maximum values of 
local stress and strain for AS-j_. An unnotched specimen was loaded 
in uniaxial tension to obtain the stress- strain curve 0A. Point A 
corresponds to the intersection of the curve for equation (10‘) and the 
stress- strain curve, and, as a result, represents- the only solution 
(Aoj_ and Ae^_) of equation (lO) consistent with the material 
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stress-strain behavior OA. Thus point A represents the estimated 
maximum values for local stress .and strain during the first quarter 
cycle. 

For the second excursion of nominal stress IS 2 point A 'was 
taken on the initial state. Correspondingly, equation (10) was plotted 
relative to point A for ASg. This curve for equation (10) lies in 
the third quadrant. The unnotched specimen was unloaded from point A 
and loaded into compression until the resulting stress-strain curve 
intersected the limiting curve for ASg. This intersection point C 
represents the estimate for the minimum values of local stress and 
strain for this cycle of nominal stress. 

Point B in figure 4(b) represents the half -cycle residual stress- 
strain state and was found by plotting equation (10) relative to A 'for 
the removal of £Sp from the specimen. Point D, representing the full- 
cycle residual condition, was found in a similar manner. 

The local stress-strain curves for cycling beyond the first cycle 
were found by repeatedly applying the procedure used to find the 
extreme points A and C. ' 

, Because this procedure was applied in a continuous cycle-by-cycle 
manner, the history dependence of local behavior was automatically 
simulated in these tests with unnotched specimens. In addition, the 
influence of the cyclic strain^hardening behavior for the material used 
in this study was inherently incorporated into the local stress-strain 

* j 

results. 

The unnotched, specimen configuration shown in figure 5*5 was used 
throughout this study. These specimens were made from 2024-T3 aluminum 
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alloy sheet, and were all oriented such that their longitudinal axes 
were parallel to the rolling direction of the sheet stock. 

The hydraulic closed-loop, servo-control testing system shown in 
figures 3*6 and 3-7 was used for cycling the unnotched specimens between 
the calculated limits for local stress and strain. While this system 
was capable of controlling either load or strain, ' neither control mode 
was used for this phase of the study. Flow control valves were 
adjusted to limit the strain rates in the specimens, and the loading 
controls were first set for capacity tensile loading. Because of a 
limited strain rate,, the stress-strain curve was slowly generated. 

When this curve intersected the limiting curve for maximum local .stress 
\ and strain (see fig. 3 - *0 the loading controls were manually reversed. 

, * t 

Similarly the loading controls were manually reversed when the stress- 
strain .curve intersected the limiting curve for minimum local conditions. 
This procedure was repeated -for each half-cycle throughout the analyti- 
cal' control test. ’Guideplates shown in figure 3* 8 were used to prevent 
buckling during the compressive loading. 

The stress-strain curves obtained in this phase of the study were 
recorded using an x-y recorder shown in figure 3«7* The load (stress) 
input to this recorder was taken from the dynamometer shown in 
figure 3*7* The strain input was provided by the strain-gage exten- 
siometers also shown in this figure. The pair of extensiometers, shown 
again in figure 3* 8, were used to determine the average strain over the 
1-inch test section of the specimen. In addition, to reduce strain 
errors due to bending, an average output from the two extensiometers 


was used. 


x— y recorder 



Figure 3-6 — Block diagram for hydraulic servo-control testing system. 
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Figure 3.7 — Photograph of apparatus used for analytical-control tests. 
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3« 3 Companion-Specimen Method 

To evaluate the accuracy of the analytical-control procedure, local 
stress and strain were determined experimentally by a companion- specimen 
method. Local strain limits are found experimentally for each loading 
cycle on the notched specimen and, these limits are used to control the 
cyclic simulation of local behavior with unnotched companion specimens. 

The notched specimens tested in this phase of the study were the 
sheet aluminum specimen containing a central hole, as shown in figure 3- 3* 
The unnotched specimen configuration shown in figure 3 « 5 was used for 
the companion specimens. 

To synchronize the notched specimen and the corresponding 
companion- specimen tests, the control system shown in figure 3*9 was 
developed. Loading on the notched specimen was controlled by a closed- 
loop servo-control system also presented in figure 3- 9* 

As the notched specimen was loaded, the averaged output from 
two small strain gages mounted at the local points and outputs from 
two similar strain gages on the companion specimen were compared at the 
servo-loop summing point. The difference between these two signals 
served as an input to a second servo-control system that controlled 
the loading on the unnotched specimen. This second control system 
maintained the strain in the companion specimens equal to that 
experienced at the local point. A dual-channel x-y recorder was used 
to simultaneously record the loads on both specimens and the common 
local strain. Photographs of the loading system and the servo controls 
are shown in figure 3*10* 
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Figure 3*9 — Block diagram of strain— coupled servo— control system. 
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Guideplates shown in this figure were used to prevent specimen 
buckling under compressive loads. 

3* k Test Program 

As previously mentioned, the objective of the local stress- strain 
phase of this study was twofold. An evaluation of the analytical- 
control technique was sought, and, in addition, 'the local stress- strain 
results- from this evaluation were to be used in the analysis of crack- 
initiation lives. Loading sequences were selected with these objectives 
in mind. 

To investigate the adequacy of the analytical- control method in 
accounting for cyclic strain hardening, two types of constant-amplitude 
loading were considered. Simple two-level loading and a random sequence 
of loading -were selected as examples of variable- amplitude loading. To 
demonstrate the use of local stress and' strain for analysis of crack 
initiation, the loading sequences and magnitudes were selected to 
produce effects on initiation behavior due to local residual stresses. 
Details of these loading conditions follow. 

3-b . 1 Constant-Amplitude Loading 

Two types of constant-amplitude loading were studied. The first 

» 1 k 

of ‘ these, shown in figure 3*ll(a), was completely reversed (R = -1) 
loading. For the notched specimen configuration and materials 

i 

t ’ 

investigated, this loading produced local inelastic behavior for both 
tensile and' compressive loading. The repeated (-R = 0) loading shewn 

V *. 

in figure ’ 3.11(b), created only elastic stresses throughout the notched 
. specimen. As a result, this loading condition was of little interest 

t 
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in this phase of the study but played an important role in the fatigue 
crack- initiation study and therefore was included here for completeness. 
5.4.2 Variable-Amplitude Loading 

The above mentioned two types of loading were combined in two-level 

t 

loading sequences shown in figure 5*12. For these conditions residual 
stresses are produced during load reversal. These residual stresses 
T affeet the fatigue resistance of the notched specimen during the second 
level of loading. 

i 

*■ 

, The fatigue behavior corresponding to these two types of loading 

. . v 

will be discussed in the next chapter. In this chapter these loading 

sequences were used to. demonstrate the role of high stresses in creating 
• 1 

local residual stresses that alter .subsequent local stress-strain 

' t 

! conditions and to demonstrate the accuracy of the analytical-control 
method for calculating these effects. The random stress sequence in 

* 9 

figure 3.15 was also investigated to illustrate the adequacy of the 
analytical- control method for arbitrary loading sequences and magnitudes. 

Local stresses and strain were found for each of the loading 
conditions in figures 3 * 11 , 3*12, and 3*13 by both the analytical- 
control and the companion- specimen methods. These results are 
presented and discussed in the following section. 

3- 5 Experimental Results 

The strain-coupled servo-control system was used to conduct 
companion- specimen tests for each loading condition. A typical 
recording from the dual channel x-y recorder is shown in figure 3*l4. 
Local stesss and strain are plotted in figure 3*l4(a) corresponding to 
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Figure 3*12 — Two-level nominal stress, sequences 
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Figure 3.14 — Typical recording from companion-specimen tests. 
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the nominal stress and local strain in figure 3-l 1 t-(b). Recordings 
■were made throughout the first 50 load cycles. Maximum, minimum, and 
residual stresses and strains were determined from these recordings 
for each type of loading. These salient local stresses and strains 
were plotted for 50 cycles of loading and -trill he discussed in the 
following sections. 

In addition, elastic stress-concentration factors were determined 
experimentally for each test specimen. Prior to the strain- coupled 
operation, strain-gage readings were taken at each of 10 nominal 
stress levels, as each notched specimen was loaded to 80 percent of 
the local yield stress. The experimental elastic stress-concentration 
factor K eX p for each .specimen was the average of the 10 values found 
for each case. These average values based on net nominal stress varied 
from 2.5^ to 2.59 with an overall average' value of 2-57 which agreed 
with the value given by Howland [I2J for this specimen configuration 

(fig. 3-3). 

3.5.I Constant-Amplitude Loading 

Maximum, minimum, and residual stresses and strains for the two 
cases of constant-amplitude loading are presented in figures 3-15 j 

\ t 

■ 1 3-3-7.; 3- 18, and 3* 19* The dashed curves also shown on these figures 

<r 

will be discussed later. 

- ; - The local stress of 51 hsi shown in figure 3-15 f° r S = 20 ksi- 

i • equals the raonotonic yield stress for 202h-T3 aluminum alloy as 
illustrated in figure 3*16. Consequently the local stress-strain 
behavior was initially elastic for this loading condition. 
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Figure 3,15 — Local stress history for constant-amplitude loading with S = 20 ksi 
. ** ; and R = 0. 
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Furthermore, the cyclic stress-strain curve, also presented in 
figure 5.16 (from [l5[]) indicated this local "behavior to he stable 
for repeated cycling, "which agrees with the results presented in 
figure 5*15* These results could have been calculated directly from, 
figure 3.16. However, the companion- specimen test was conducted in 
the interest of completeness and also to assess the performance of' the 
strain-coupled servo-control system. In addition, this local stress 
history for S = 20 ksi will be used as a reference in subsequent 
sections to describe the local effects of prior high loading on 
subsequent load cycling at S = 20 ksi. The local effects will in 
■ turn be', used to explain the crack- initiation behavior during fatigue 
tests of the notched specimens. 

The application of S = kO ksi produced local reversed plasticity 
during each cycle as shown in figure 3* 17* The cyclic strain-hardening 
behavior shown in figure 3.l6 for 202^~T3 aluminum alloy was also 
evident in figure, 3* 17* The ma ximu m and minimum local stresses 
increased initially with cycling but stabilized in about 10 cycles. 
Residual stresses decreased in magnitude corresponding to the increases 
in maximum and minimum stresses. 

The strain history in figure 3-18 again illustrates that the local 
behavior for S = 20 ksi was "both elastic and stabilized. The" strain 
range in figure 3.19 for S = ItO ksi, however, decreased rapidly for 
the first few cycles before approaching, a stabilized condition in 
about 10 cycles. 



Figure 3. .17 — Local stress history for constant-amplitude; loading with S = 40 ksi 
and R = —1. 
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Figure 3 .I 8 — Local strain history for constant-amplitude loading with S = 20 ksi 
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Figure 3.19 - Local strain history for constant-amplitude loading with S = 40 ksi 
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3*5*2 Variable Amplitude Loading 

As previously mentioned in the discussion of the test program, the 
two constant-amplitude conditions were combined to form two distinct types 
of two-level loading. The high loading (s = Uo ksi) was applied first and 
was followed by load cycling at S = 20 ksi. The local plasticity pro- 
duced by the high load cycling resulted in the creation of residual 
stresses that altered the local stress- strain conditions for subsequent 
low load cycling. The high loading was applied for 10 cycles and ended 
on both tensile and compressive half-cycles to produce, respectively, 
t beneficial (compressive) and detrimental (tensile) residual stresses. 

For the case -of beneficial two-level loading a compressive residual 
stress of ‘30 ksi' was left at the end of the high loading phase, 
figure 3*20. This residual stress acted as the minimum local stress 

’ ^ i 

during the second phase of loading, and decreased slightly during the 

- J , ; ' . . 

/ t , , ’ t 

first hQ cycles '.of loading at the low stress level. The maximum local 
stress corresponding to the loading at S = 20 ksi was initially 22 ksi, 
•but also decreased slightly as the cycling was continued. This decrease 
in maximum and minimum stress shown in figure 3*20 was believed to -be the 
result of strain-gage zero shift, and thought not to be representative 
of actual local stress behavior. This will be further discussed later 
in this report. 

A tensile residual' stress is shown in figure 3*21 for high loading 
ending in compression. For this case the minimum local stress was 7 ksi 
and the maximum local stress was 58 ksi for the cyclic loading at 
S = 20 ksi. Again for this case, a slight decrease occurred in extreme 
stresses as the cycling continued. 


local stress, a , ksi 



Figure 3.20 — Local stress history for beneficial prior loading. 
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Figure 3.21 — Local stress history for detrimental prior loading, 
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The local strain histories for the two types of two-level loading 
are presented in figures 5*22 and 3*23* The local strain ranges 
decreased during the 10 cycles of high loading, but remained virtually 
unchanged for load cycling at the low stress level. The tensile 
residual strain at the end of 10 cycles in figure 5*22 produced the 
compressive residual stress shown in figure 5*20. Correspondingly, the 
compressive residual strain from figure 5*23 caused the tensile 
residual stress in figure 3*21. 

' , Maximum and minimum local stresses and strains for the random 
stress sequence in figure 3*13 were determined by the companion- 
specimen method and are presented in figure 3*24., The notched specimen 
buckled under the first application of -50 ksi, in spite of the guide 

plates. For this reason the stress and strain histories were- terminated 

1 " * 

after six cycles. Nevertheless, these experimental results are pre- 
sented and will be used in the next section to assess the accuracy of 

f . 

the analytical- control method for random loadings. 

3. 6 Comparison of Calculated and Experimental Results 

For an evaluation of the analytical-control method, this method 
was applied for each previously discussed loading condition to calculate 
the local stress-strain behavior from unnotched specimens. A typical 
recording of local stress- strain curves obtained by this approach are 
shown in figure 3*25. The calculations, of local behavior by analytical 
control were continued for 50 cycles as in the companion- specimen tests. 
The calculated results are presented as local maximum,, minimum, and 
residual stresses and strains for comparison with similar experimental 
data., The analytical-control results are represented by the dashed 
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Figure 3.23 — Local strain history for detrimental prior loading. 
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curves plotted together with the companion- specimen results in 
figures 5*15 and 3 * 17 ' through 3-23* 

5.6.1 . Constant -Amplitude Loading 

The ' calculated local stress history for S = 20 ksi was entirely 
elastic., as 1 expected, and agreed very closely with the corresponding 
experimental results in figure 3.15* This case, of course, did not 
represent a critical evaluation of the analytical-control method, but 

was studied to obtain a reference for local conditions for an evalua- 

( 

tion .of loading sequence effects during two-level loading. The constant- 

amplitude loading at S = 40 ksi provided a better test of the 

analytical-control approach. The calculated- results for this case 

also correlate very well with data in figure 3.17* This correlation 

demonstrates the ability of the analytical- control .method to account. 

for local cyclic plasticity and to inherently incorporate the cyclic 

strain-hardening behavior in calculations of local cyclic stress. 

As expected, the calculated and experimental strain histories 

for S = 20 ksi agreed closely, as illustrated in figure 3*l8. In 

contrast, the calculated local strains in figure 3-19 for S = 40 ksi 

were somewhat in error. 

\ 

3*6.2 Variable-Amplitude Loading 

Local stress histories calculated for both types of two-level 
loading correlated remarkably well with the experimental results in 
figures 5-20 and 3*21. The stabilized calculated results for S = 20 ksi 
support the assumption that the slight decreases in .local experimental 
stresses were due to small strain-gage errors. As in the case of 
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constant - amplitude loading^ the calculated local strains were noticeably 
larger than the experimental values in figures 3.22 and 3.23. The 
maximum and minimum stresses and strains for the random loading sequence 
are presented in tables I and II. The corresponding experimental 
results from figure 3«24 are also listed in these tables and agree very 
well with the calculated results. 


TABLE I.- COMPARISON OF CALCULATED AND EXPERIMENTAL 
LOCAL STRESSES FOR RANDOM NOMINAL STRESS SEQUENCE 


Cycle 

Number 

Maximum local stress, 
ksi 

Minimum local stress, 

ksi j 

■ Analytical- 
control 
method 

Companion- 

specimen 

method 

Analytical- 

control 

method 

Companion- 

specimen 

method 

1 

54.1 

53-8 

' -55.3 

-54.1 

2 

50.8 

50.0 

0.0 

-2.0 

.3 

56.1 

55.8. 

-19.4 

-19.8 . 

4 

58.I 

58.6 

-37.1 

-34.8 

5 

• 38.7 

41.2 

t-64. 6 

-63*9 

6 

-11.3 

-10.9 

-37.5 

-37*7 

7 

48.8 

48.3 

-67.0 

(specimen 

8 

55.6 

♦ 

-67.0 

buckled) 

9 

11. 7 


-67.0 


10 

21.0 


-67.0 


ll 

'65-O 


-10.5 


12 

15.3 


-69.O 


13 

58..I 


-41.2 


14 

63.0 


-64.6 

1 t 

15 

67.O 


-50.0 


16 

50.0 


-75.1 


17 

69.4 


-77.1 


18 

58.1 


-64.6 


19 

67.O 


-75*9 


20 

• 47.2 


421.0 


21 

70.2 


-61.3 


22 

70.2 


-72.6 


23 

69.4 


-73*0 


24 

29.9 


' 4-4.0 


25 

71.0 


-46.0 


2 6 

29.9 


-67.0 


27 

70*2 


-67.8 


28 

52.5 


4-26.6 


29 

■ 70.2 


-71. 8 


30 

48.4 


-0.8 

■ 
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TABLE II.- COMPARISON OF CALCULATED AND EXPERIMENTAL 
LOCAL STRAINS FOR RANDOM NOMINAL STRESS SEQUENCE 


Cycle 

Number 

Maximum local 
stress 

Minimum local 
stress 

Analytical- 

control 

method 

Companion- 
specimen ' 
method 

Analytical- 

control 

method 

Companion- 

specimen. 

method 

1 

1.07 

0.89 

-I.06 

-1.01 

2 

0.14-7 

0.40 

- 0.03 • 

-0.11 

3 

1.00 

O.89 

0. 25 

0.13 

4 

1-57 

1-75 

0.49 

O.65 

5 

1 . 2 k 

1.42 

-1.80 

-1.66 

6 

- 1.33 

1.19 

-I.56 

- 1.42 

7 

-0. 38 

-0.36 

- 2.37 

■ ( specimen 

8 

-0.50 


- 2.37 

buckled) 

9 

- 1 . 6 k 


- 2.37 


10 

- 1-55 


- 2.37 


li 

1.514- 


0.75 


12 

1.00 


■ - 1.71 


13 

0.11 


-O.95 


14 

0.60 


- 1.23 


15 

1.92 


0.54 

- 

1 6 

1.57 


- 2.55 


• ’ 17 

1. 914- 


-2.50 


18 

-0.81 


-2.11 


•19 

0.29 


-2.60 


20 

-1.31 


- 1.55 


21 

1.67 


-0.10 


22 

. 1.67 


- 1.23 

« 

23 

1.69 


-1.22 


2 k 

. -0.20 


- 0.45 


25 ' ■ 

2.98 


-1.61 


26 

2.37 


0.70 


27 

2.99 


0.73 


28 

2.07 


1.82 


29 

3.13 


0.28 


’ • ^30 
• ' * 

1.57 


i.o 4 

- 












4. CRACK INITIATION AT STRESS CONCENTRATIONS 


if-. 1 General 

The analysis of fatigue in service is usually extremely complex as 
a result of the very localized nature of fatigue. Fatigue cracks 
initiate at the most severely stressed area of a machine part or 
structure and this area of high stress concentration is characteris- 
tically very small. A prediction of fatigue behavior must be based 
upon the fatigue behavior of this highly stressed local zone. In 
addition, fatigue predictions must account for the growth behavior of 
cracks that originate and grow from this local zone, but in this chapter 
attention, is focused on crack initiation. 

In principle, fatigue crack initiation at a stress concentration 
can be predicted by calculating local stresses for each cycle of loading 
and by estimating the corresponding " local" fatigue damage from a 
cumulative damage theory and a fatigue crack- initiation diagram for 
unnotehed specimens. This approach is complicated by local plasticity 
and attendant residual stresses, by interactions of one local stress 
level with the damage caused at another local stress level, and by a 
size effect. 

Local plasticity and the resulting residual stress behavior were 
illustrated in figures 3.20 and 3*21 and discussed in the previous 
chapter. From a fatigue point of view, residual stresses influence 
local damage-accumulation rates by altering local mean stresses. 

Because fatigue damage increases with mean stress, the prior loading 
in figure 3* 21 should have a detrimental influence on the local fatigue 
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behavior under S = 20 ksi loading. On the other hand, the prior 
loading for figure 3*20 should improve the fatigue performance under 
the low level. These nominal stress-interaction effects may be accounted 
for only by accurate calculations for local stresses. As indicated 
in Chapter 2, such calculations are' quite complicated and, as a result, 
approximate methods are usu ally required. 

In addition to the stress interaction due to residual stresses, a 
second type of stress interaction is produced by damage acceleration. 
Dolan and Corten Qtjj] and Manson, Freche, and Ensign [W] described 
damage acceleration in terms of microscopic crack initiation caused by 
a high stress level that enabled low stress cycles to readily propagate 
the initiated microcrack. A considerable fraction of the life would be 
required in constant-amplitude, cycling at the low stress level to 
create a similar microcrack. Consistent with this explanation, a few 
high stress cycles, have been observed to noticeably accelerate the 
damag.e accumulation for subsequent cycling at a lower stress level. 
Freudenthal and Heller described damage acceleration in terms of 

the increased crack-propagation rate at each stress level caused by 
interaction with all higher stress levels. In both explanations 
damage acceleration was related to microscopic plastic deformation at 
the high stress that enabled the lower stresses to do more damage th an 
they could otherwise do. 

Both types of stress interaction were simulated in unnotched 
specimens by the analytical-control procedure from Chapter 3 - By 
continuing an analytical-control test until fatigue cracks initiate. 
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these interaction effects can he automatically incorporated into the 
resulting estimate for fatigue crack initiation. Furthermore, the 
fatigue damage accumulation in the unnotched specimen simulates that 
for the material at the notch root, eliminating the need for a cumula- 
tive fatigue damage theory. Estimates of crack-initiation periods were 
found by this analytical- control approach and are compared with data 
in this chapter. 

The analytical-control procedure for estimating crack initiation 
does not account for the size effect that causes Kf < Kj. This size 
effect results from the fact that when local stress is simulated in 
unnotched specimens, a large volume of material is highly stressed, 
although only a small volume is similarly stressed at the actual stress- 
concentration site. From a statistical point of view, the large volume 
of material in the unnotched specimens contains a larger microscopic 
flaw than the smaller volume .at the stress-concentration site. This 
larger flaw in the unnotched specimen grows to a macroscopic size faster 
than the smaller flaw at the stress concentration and, as a result, 
for the same stress history, the crack- initiation period for the 
unnotched specimen is shorter than for the stress concentration. 

Wetzel [24] replaced Kf hy Kf in 'the Weuber equation '(9) to account 
for this size effect and consequently, to improve fatigue estimates. 

i 

Although for cases of constant -amplitude, Wetzel has improved estimates 
of crack initiation. < This ’ approach leads to large discrepancies 

* i , 

between simulated and actual local stresses,’ especially for high 

1 * 

stress-concentration factors where Kf is usually significantly smaller 


t 


{ 


than Rj. These discrepancies in local stress would alter both of the 
.previously explained stress interactions for variable-amplitude loading. 

t ' * > 

*r i * « . 

For this reason no correction was made .for'' size effect in the present 

f 

study. 

4.2 Observed Crack- Initiation Periods 

The sheet specimens containing a central hole from Chapter 3 were 
tested under fatigue loading until fatigue cracks developed at the 
stress-concentration site. Tests were conducted for the two levels of 
const ant -amplitude loading and "the two combinations of two-level loading 
considered in Chapter 3. 

Microscopic cracks usually develop at an early stage of fatigue 
life. However, a considerable fraction of the total fatigue is often 
required to propagate these small cracks to a macroscopic size. From an 
engineering point of view, very small cracks (that do not alter the load 
carrying capacity of a structural member) may be ignored. From this 
point of view, the crack-initiation life is often defined as the number 
of cycles required to produce a macroscopic crack, detectable by con- 
ventional inspection methods. For the present study, crack- initiation 
life was defined as the number of cycles necessary to create a 
0.03-inch crack at the stress -concentration site. This "initial” crack 
size of 0.03 inch was selected from crack propagation considerations 
(see section 3*2)-, but nevertheless, is representative of the minimum 
size detectable by inspection methods currently In service. 

Because of difficulties in direct detection of small cracks 


(0.03 inch) at the stress-concentration site for each fatigue test, an 


i 
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indirect approach was used to establish crack- initiation lives. The 
total fatigue life, Nf, of notched specimens was assumed to consist of 
two parts: the number of cycles required to produce the 0.03-inch crack, 

SF 0 , and the number of cycles necessary to propagate this crack to 
failure, N C p. As part of the crack propagation study in Chapter 5, 

N C p was experimentally determined for the notched specimen configura- 
tion for each loading sequence used in the present study of crack 
initiation. For each fatigue test in this chapter, the crack- initiation 
life, N 0 , was found by subtracting the corresponding constant value- of 
w c p from the observed total life, These N cp values were: 

*K)00 cycles for the case of no prior loading. (S = 20 ksi'), 34,000 cycles 
for beneficial prior loading, and 3000 cycles for detrimental prior 
loading.. Errors'- associated with the assumption of a constant N _ • 

value for each ‘test condition were believed to be small, since for all 

1 < < * 

test conditions these N cp values- were less, than 5*6 percent of the 

* * * 
corresponding geometric, mean'" values of. 

* ' t 

Total fatigue lives and the corresponding crack- initiation periods 
are presented in table III -and figure 4.1 for each of the four test 
conditions. In all tests cracks are initiated at the zone of highest 
stress concentration. In addition to 12-inch-wide specimens, 4-inch- 
wide specimens were used as indicated in table III and illustrated in 
figure 4.1 by open symbols. The constant-amplitude loading conditions 
resulted in initiation lives with geometric means of 331 cycles for 
S = 40' ksi and 129,000 cycles for S = 20 ksi. The high- low combina- 
tions of these two stress levels produced different fatigue lives. 


TABLE III. - FATIGUE CRACK- INITIATION DATA FOR BOTCHED SPECIMENS 


Botched 

Test condition specimen 
* number( a ) 


S = 40 ksi 
R = -1 


(B cp = 12 cycles ) 



Beneficial 
prior loading 

(N C p = 34,000 cycles) 


, Detrimental 
prior' loading 

(N C p = 3000, cycles) 


108-1 
98 • 

-1TP-P 


Geometric mean 
of crack initiation 
lives, N 0 , cycles 


97-1 

109 

108-4 

112-1 

93-5 



93 

* 103' ’ 
l68 ; 
112 . 
93-1 
, 103-1 
103-3 


100 267, 300 

102 • 303,100 

93-3 ’ 590,.200 

112-3 910,000 

93-4 1,003,100 

97-4 1,063,000 

108-2 1,071,000 


54,200 

65.000 

66.000 

77,300 

78. 000 


233,300 

269.100 
556,200 
876,000 

969.100 

1.029.000 

1.037.000 


51,200 

62.000 

63.000 
74,300 

75.000 



^'Hyphenated numbers indicate 4- inch-wide specimens with 0.672-inch- 
diameter' hole . 

(b)Crack- initiation lives found by subtracting, crack propagation period 
N cp (see first column) from total life Bp. 
























< 


(No- vertical scale) Material: 202MT3 #—12 in-wide specimens 



Figure k.l — Observed crack initiation behavior for notched specimens. 
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As expected from the analysis of local stress-strain conditions, the 
loading sequence in figure 3.12(a) resulted in a fatigue life larger 
than the reference case of constant-amplitude loading, while the 
sequence in figure 3» 12(b) produced a shorter life. That is, the prior 
loading in figure 3* 12(a) was beneficial and the prior loading in 
figure 3* 12(h) was detrimental. To emphasize this stress interaction 
the linear cumulative damage rule (ref. 29) 



(18) 


was applied where M is' the number of stress levels, np is the 
number of applied cycles for the ith stress level, and N-j_ is the 
constant-amplitude fatigue life 'for that stress level. For the case 
of beneficial prior loading 


V _ fl + ^2 _ _10_ + 610,000 
L % Up U 2 331 129,000 


4.77 


and for detrimental prior loading 


^ ^ = _1Q_ + 63.U0Q _ 0 =* 

% n 2 331 129, 000 * 

These results illustrate the nonlinear damage accumulation caused 
by residual stresses not accounted for by the linear cumulative damage 
rule. The difference in the last half -cycle of prior loading for these 
two cases (compare figs. 3* 12(a) and 3- 12(b)) caused approximately an 
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order of magnitude difference in crack- initiation lives. For service 
loading sequences, which in many cases are similar to the random 
sequence of loads in figure 3-13j residual stresses are being created 
and eliminated throughout the entire crack- initiation life. Fatigue 
predictions for such cases must he based on the analysis of local, 
stress- strain conditions to account for the complex residual stress 
effects. The application of the linear cumulative damage rule based 
on nominal stresses leads to meaningless results for these service 
loading conditions. 

^.3 Prediction of Fatigue Crack Initiation 

To obtain estimates of crack- initiation periods for the notched 
specimen, the analytical-control procedure presented in Chapter 3 was 
continued until the unnotched specimens failed in fatigue. As 
previously illustrated (figs. 3»17 j 3*20, and 3*21), local stress-strain 
conditions stabilize ‘for repeated load cycling. For convenience, 

* ‘ f 5 

after stabilization occurred the testing procedure was changed from 
analytical control to stress control and continued in the stress 

* * ’ t * 

control mode until failure '-occurred. Since crack- initiation life for 

* i * * 

the notched specimens was, defined as the number of cycles required to 

t 

produce a 0.03-inch crack, these tests with unnotched specimens should 

* i 

‘ * 

have been continued until a crack 'of this same size developed. However, 
preliminary tests showed that these unnotched specimens failed statically 
for small crack lengths, comparable to the 0.03- inch "initial" size 
crack. Consequently,, for convenience crack propagation was not 
monitored in these tests, and the estimates for crack- initiation lives 
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were taken as the number of cycles for complete failure of the 
unnotched specimens. 

The errors introduced in N 0 predictions by these simplifying 
assumptions were believed to be small compared to the effects of load 
interaction. Both the residual stress and the damage acceleration 
stress interactions were accounted for by the • analytical-control 
fatigue tests. It is believed only the unaccounted for size effect 
resulted in significant errors. 

Analytical-control tests were conducted for the four test condi- 
tion's and resulting estimates’ for N 0 are presented in table IV 
and figure 4.2. The beneficial and detrimental effects of residual 
stress are reflected by the results' in figure 4.2. These results are 
further discussed in' the next section and are compared with the observed 
crack- initiation periods. 

4. 4 Comparison of Predicted and Observed Results 

For an assessment of the accuracy of the analytical-control method 
for estimating fatigue, the observed and predicted results from 
figures 4.1 and 4.2 were combined in figure 4.3. The predicted crack- 
initiation periods reflect the general trends displayed by the observed 
results, however, the previously mentioned size effect and stress 
interactions combine to produce conservative predictions. These effects 
are each briefly discussed in the following analysis of the discrepancies 
between predicted and observed crack- initiation periods. 

The analytical-control method as presented in this report does not 
account for the size effect on the fatigue crack- initiation behavior. 
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TABLE IV,. - PREDICTED CRACK - Hi IT IATION PERIODS FROM 
ANALyTICAL-CONTROL TESTS WITH TJNNOTCHED SPECIMENS 

t ! , i 


Test condition 
for notched 
specimen • 

Unnotched 

specimen 

number 

» 

Fatigue life 
. from analytical- 
control tests, 
cycles 

Geometric mean, 
of fatigue lives, 
cycles. 

■ 

1 

i 

96 


1 ' 


166 

155 

S = 40 ksi 


232 


R = -1 

■ 




• G53 ' -• 

• 19, 200 



G51 

22,700 


S = 20 ksi 

N53 

43,500 

35,100 

R = 0 

112-A 

48,100 



L51 

58 , to) 

1 


M51 

140,000 



A53 

181,000 

191,000 

Beneficial 

F53 

273,300 


prior loading 





C51 

17,300 


Detrimental 

' D 51 

23,100 

22,400 

prior loading 

F51 

28,800 


















(No vertical scale) 
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E = -1 
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E = 0 
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prior loading 
(see fig. 3*124) 
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Figure 4.2 — Predicted fatigue crack initiation periods for notched specimens 
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prior loading 
(see fig. 3.12B) 
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Figure 4 . 3 — Comparison of 
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The unnotched specimens* used to estimate W 0 contained a much larger 
volume of highly stressed material than notched specimens with the 
result that observable fatigue cracks initiated earlier in the 
unnotched specimens, as previously explained. This effect is attenuated 
somewhat by local plasticity because a larger volume of material 
(within the local plastic zone) experiences the high local stress. 

This trend was present in the predicted W Q values in figure 4, 3 for 
constant-amplitude loading. The estimate of N 0 for S = 20 ksi was 
a jRP r °ximately one-fourth the corresponding observed value (see 
tables III and TV for .actual values.). In contrast, the estimate for 
S = ko ksi, which produced local plasticity, was nearly one-half the 
observed N 0 . To illustrate the effect of plasticity on the size of 
the highly stressed local zone, the elastic stress distribution in 
figure 3.2 was used in equation (ll) together with the stress-strain 
curve in figure 3.16 to obtain the stress distributions in figure 4, 4. 
The stress distributions in figure 4.4 for prior loading show that the 
size effect should also be diminished for both cases of two-level 
loading. For these two test conditions the predicted I'T 0 values were 
approximately one-third as large as the observed N 0 values (see 
fig. 4.3 and tables III and XV). 

The size effect is usually expressed in terms of Kg>, the ratio of 
fatigue stresses for- unnotched -and notched specimens at a given fatigue 

I *, ^ , j * * 

life. Auxiliary fatigue tests were conducted- with unnotched specimens 

to establish the value of Kp for this study. Results from these tests 

t ' > ' , V 

are listed in table V and are plotted in, figure 4. 5 to create a segment 
• * ‘ * • * 
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of the S-N curve for unnotched specimens. On this curve a stress of 
4j5*5 ksi corresponds to the geometric mean of the lives for notched 
specimens at S = 20 ksi. The ratio of stresses at this life yields 
% = P.l8 in contrast to Kp = 2.57 from Chapter 3 . This value 
of Kf could have been used in equation (10) in place of % to 
improve estimates of h 0 for constant-amplitude tests as suggested. by- 
Wetzel [24], however, this substitution would have introduced noticeable 
errors in the local stress- strain calculations for the cases of two- 
level loading. These errors probably would not have seriously affected 
the N q estimates for the present study because the above Kp is 
85 percent of Kp. However, for high Kp values, the difference 

i 

between Kp and Kp is usually quite large and Wetzel' s approach 

would be unacceptable for variable- amplitude loading, because of 

inappropriate stress-interaction effects., A general correction 

*+ 

procedure for the size effect must be accomplished through appropriate 
* 1 * 1 

adjustments of the , estimated N 0 found by 'cycling unnotched specimens 

* , ! 

through the, most accurate local stress (or strain) sequence obtainable. 

’ • - ’ ■ • ’ 

To enable an evaluation of the effects of stress ' interactions on 

the analytical- control tests, additional constant-amplitude- test results 

in table y and figure 4.6 were obtained. These tests were designed 

to illustrate separately the residual stress effect and the damage 

acceleration effect. Unnotched- specimens were cycled between the 

constant-amplitude local stress limits corresponding to S = 20 ksi 

for the two-level tests to directly assess the effects of residual 

stress. These stress limits were determined from analytical-control 



□ — Constant— amplitude tests 
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Figure k.6 — Effects of mean stress and prior loading on fatigue of unnotched specimens 
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results in figures 3*20 and 3*21. The analytical- control estimates, 
for W 0 from table ( IV were also presented in figure 4. 6 for compari- 
son. These analytical-control results reflect the effects of both 
types of stress interaction. Consequently, the difference between 

these results and the constant-amplitude results is a measure- of the 

, * 

i ' * 

damage acceleration' effect for each test condition. For comparison the 
analytical-control results corresponding to constant-amplitude S = 20 ksi 
are also presented iri figure 4.6. 

The residual stress effect increased the geometric mean of the 
estimate for N 0 from 33,100 to 302,000 cycles in the case of bene- 
ficial prior loading. However, damage acceleration caused by the 
10 cycles of high load reduced this estimate to 191,000 cycles as shown 
in figure 4.6. The net effect was an increase in life due to the 
decrease in local mean stress caused by the 10 cycles of prior loading. 
These effects agree generally with results obtained by Topper and 
Sandor in a parametric study of the effects of mean stress and 

prestrain on fatigue damage accumulation in 2024-T4 al uminu m alloy. 

For the case of detrimental prior loading, corresponding constant- 
amplitude tests produced a mean life of 31,600 cycles and simulation 
tests resulted in a mean life of 22,400 cycles. The residual stress 
effect decreased the life estimated from 35,100 to 31,600 cycles, and 
damage acceleration further decreased the simulation results to a mean 
life of 22, 400 cycles. These effects are also in general agreement with 
results' obtained by Topper and Sandor '[46] and, in addition, further 
illustrate the relative magnitudes of the two stress interactions. 
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In summary, the role of stress concentrations in stress-interaction 
effects has been demonstrated in this chapter by crack-initiation tests 
with notched specimens. The stress interactions considered in this 
study were shown to be .the combined result of residual stress effects 
and damage acceleration caused by large excursions of local stress and 
strain. The procedure described for predicting crack- initiation lives 
accounted for both of these effects and produced estimates that 
reflected the general stress -inter action effects. However, these 
predictions were consistently shorter than experimentally determined 
crack- initiation lives for notched specimens. This trend was attributed 

T f 

to the large difference in volumes of highly stressed material for the 
notched specimens and the unnotched specimens used to obtain the 
predictions. .Additional research is required to assess this size effect 
and to develop rational procedures for accurately predicting crack 
initiation at stress-concentration sites. 



5»1 General 


5*,' CRACK PROPAGATION HEAR STRESS CONCENTRATIONS 

* r 



The role of a typical stress, concentration in' the fatigue crack- 

initiation process was introduced and discussed in Chapter’ 4. Local 

* 

stresses at a stress-concentration site were calculated and used 
together with unnotched specimens to predict fatigue crack- initiation 
periods. To complete the fatigue analysis for this typical stress 
concentration, the behavior of these fatigue cracks was investigated 
as the cyclic loading was continued beyond that required for crack 
initiation. The initial stage of growth for a crack emanating from a 
stress -concentration site is influenced by the stress distribution 
caused by the stress concentration. However, for repeated loading the 
crack eventually attains a length for which the influence of the 
stresses due to the original stress concentration may be ignored, and 
the crack then may be assumed to propagate solely under the - influence 
of the applied nominal stress. This phase of crack growth has received 
much attention in the literature (for example (V/j , [48], and [ 49 ] ) and 
was not specifically investigated in this study. Emphasis was placed 
on crack growth under the influence of the stress concentration. This 
initial phase of crack growth from a stress -concentration site has 
received very limited coverage in the literature, as indicated in 
Chapter 2. 

The present study was undertaken to provide an experimental and 
theoretical investigation of the initial stage of crack growth from a 
typical stress concentration. For this purpose, the notched specimen 
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Figure 5.2 - Crack growth curves for constantr-amplitude loading. 
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and the fatigue loading discussed in Chapter 3 (figs. 3* 3* 3»11> and 

* 1 \ 

3-12) were used for convenience and, as a result, the study of crack 

> i * 

growth was a logical extension of the fatigue analysis presented in 

, * -v 

t t 

Chapter 4. The case of constant-amplitude loading for S, = 20 ksi 

i * ' 

• . ' . 1 

was studied to demonstrate the effect of stress concentrations on crack 

i * , v 

propagation for- completely elastic 1 local conditions and, in addition, 

•» * 

this' case served as - a reference for this evaluation of the stress 
interaction effects on initial crack growth ‘under two-level loading. 

4 * 

' » t * 

The theoretical phase of this study resulted in the development 
of a method for calculating initial crack growth for an arbitrary 
elastic stress distribution near the stress concentration. This method 
was applied for the reference case of constant amplitude, S = 20 ksi 
loading, and the calculated results were compared with observed crack- 
growth behavior. As a more critical .evaluation of the method,, calculated 
and experimental results were compared for beneficial two-level loading. 
These calculated and observed crack-growth results' are presented and 
discussed in the following sections. 

5*2 Observed Crack Propagation Behavior' Near a Circular Hole 

A specimen with a central hole containing cracks is presented in 
figure 5*1 to illustrate typical crack location and to indicate the 
two schemes for defining crack length. Crack propagation tests were 
conducted by continuing the crack initiation tests from Chapter k 
beyond the initiation phase. Crack lengths were measured during these 
tests with a machinist scale and a low magnification (lOX) microscope. 

In most tests, cracks did not simultaneously initiate on both sides of 
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the hole as shown in figure 5.1. To produce the desired symmetry, after 

one crack initiated a second crack was simulated on the other side of 

the hole ty a small saw cut with a sharp nick at its root. A Jeweler' s 

saw and a sharp razor blade were used to produce this artificial crack. 

In all cases a> fatigue-crack promptly initiated at the root of the saw 

cut and symmetrical crack' growth followed. 

* 

The crack-growth curves, including the initiation periods, are 

, X * 

’ - - % 

presented in figure 5*2 for the two constants amplitude test conditions. 

, i - . . ' * - <. * " 

For simplicity, a typical curve was plotted for each test condition. 

’ * *<■ ^ 

Local cyclic plasticity accompanied' the initial crack growth for 

i t * 

S' = ifO ksi loading, precluding an analysis of this case. In contrast, 
local behavior at the hole whs' completely elastic during the cycling 
at S = 20 ksi, and enabled use of Fracture Mechanics- theory. The 
crack-growth curve for S = 20 ksi was used as a reference for 
establishing the effects of two-level loading on initial crack 
propagation. For this purpose this curve will be subsequently replotted 
with an expanded cycles scale. 

Prior loading was demonstrated in Chapter k to be either detrimental 
or beneficial to crack initiation, depending on the sense of the 
residual stress at the end of the prior loading. This behavior i-s 
again displayed by the typical crack-growth curves in figure 5 . 5 ' for- 
two-level loading. To separate these curves into crack- initiation and 
propagation stages, the initiation period was assumed to equal the 
number of cycles required to produce a crack length of 0.05 inch 
(the smallest crack length for which crack propagation data were 


normalized crack length, ~Ljp 


1 



Number of cycles, N 




obtained for all three loading conditions). The crack lengths beyond 

this initial value were plotted in figure 5*4 against the number of 

cycles for crack propagation Ii C p. Starting at the same initial length, 

the beneficial prior loading decreased the initial crack-growth rates as 

compared to the reference case of constant-amplitude loading. The 

detrimental prior loading caused the crack-growth rates to increase. 

Beyond a crack length of approximately 0.7 inch, crack growth was 

nearly identical for all three cases, indicating the localized nature 

of the influence of the stress concentration. A procedure is presented 

in the next section for evaluating this influence. 

5, 5 Procedure For Calculating Crack Propagation Behavior Near a 
Circular Hole 

Paris, Gomez, and Anderson [49] postulated the existence of a 
unique relationship between stress-intensity factor and rate of crack 
propagation. Thus crack growth could be predicted from calculated 
stress-intensity factors and "master" crack propagation rate curves, 
characterizing the crack-growth properties of the material. Stress- 
intensity factors can be calculated from elastic analyses of bodies 
containing cracks. The "master" crack-growth rate curve can be 
established by measuring crack-growth rates corresponding to known 
stress-intensity values using simple laboratory .specimens. This general 
approach was followed by Figge and Newman [50] to predict crack growth 
for complex specimen, configurations and loadings. A generally close 
correlation was found between predicted and experimental results. 

In the present study, this general stress-intensity .approach 
was used to predict crack-growth curves corresponding to the 
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constant- amplitude loading for S = 20 ksi and the case of beneficial 
prior loading. 

5* 3-1 Calculation of Stress -Intensity Factor For Crack Emanating 
From a Circular Hole. 

Bowie Qxfj developed polynomial mapping functions for use with the 
complex stress function technique to solve the problem of cracks 
emanating from a circular hole in an infinite sheet. For the case 
of two cracks emanating from a circular hole (see fig. 5.1) Paris' and 
Sih [52] presented Bowie's solution in the form 

k = °g \/L F^j (19) 

where 0g is the gross nominal stress and L is crack length. This 
stress intensity is plotted in figure 5.5. The stress-intensity 
solution 

, k = cfg \/a = a g \/L + p (20) 

> w 

» « 

for a crack of length 2a/ without the circular hole, is also plotted 
in figure 5.5 for comparison. Beyond L/p = 0.25, the curves for 
equations (l-9) and (20) are in. rather close agreement. Consequently 
for L/p > O.25 the presence of the hole may be ignored and the 
crack-growth behavior approximated by that of a simple crack with 
a = p + L, again illustrating the localized effect of stress 
concentrations. 

Bowie's solution could have been used to predict the crack growth 
for the reference case of S = 20 ksi constant-amplitude loading, 
because the local conditions at the hole were elastic. However, for 


Stress intensity factor, K, ksi 
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the two-level loading sequence, residual stresses altered the stress 

distribution around the circular hole, and consequently changed the 

stress-intensity factors for the cracks growing from the hole, as 

evidenced by the crack-growth curves in figure 5. 4 . To account for 

the effects of residual stresses on crack growth from circular holes, 

a general stress -intensity solution was developed for an arbitrary 

stress distribution cr(x) along the plane of the crack. This stress- 

intensity analysis Is presented in the following development and is 

specialized in the next section for the elastic reference case and 

for the case of beneficial prior loading. 

' « 

i 

The analysis, of stress’ intensity for the cracked specimen and 
loading in figure 5. 6(a) was based on the separation of this problem 

1 

into those given in figure 5. 6(b) and fig. 5.6(c). The stress 
1 

distribution in figure 5* 6(b) along the plane of the crack was 
designated a(x). The’ superposition of this internal stress o(x) 

. i 

from figure 5* 6(h) and ; the -a(x) shown in figure 5*6(0), acting on 
the crack boundary, resulted in the stress-free crack surface shown in 
figure 5 * 6 (a). Because the specimen in figure 5 * 6 (b) did not contain 
a crack, the stress-intensity factor for figure 5.6(a) was equal to 
that for figure 5.6(c). The stress-intensity factor for figure 5.6(c) 
was found by an additional superposition scheme used together with a 
boundary collocation technique. Irwin ' 53 presented the Westergaard 
stress function 


Z x (z) = 


2 Pz 

fa 2 - b 2 ' 

jtttz 2 - b 2 ) 

z 2 - a 2 


1 - 1/2 


( 21 ) 


for .the problem shown ip. figure 5.7 for symmetrical concentrated loads 

on a crack surface. 'This- solution was used to represent the cr(x) 

from figure 5*6(c) by a statically equivalent system of forces shown 
! * 4 

typically in figure 5.;'8. -The radial stresses cr rr and shear stress 
*T r 0 along., the circle of unit radius were found for each concentrated 
force from the following equations, developed in Appendix A. 


cr rr = jjleZ^ - yI m ZjJcos 2 e + jseZ^ + y^ n ZjJsin 2 0 
- 2yReZ£ sin 8 cos 8' 


and 


> ( 22 ) 


T r0 = yj^m 2 ! sin 0 cos 0 “ ReZ i( c °s 2 e - sin 2 0)TJ 

J 

The c rr and r r0 stress components for all concentrated forces were 
summed to obtain the stress distributions along the unit circle, 
corresponding to the cr(x) in figure 5.6(c). Following a boundary 
collocation procedure presented by Newman a second set of 

and T r 0 was expressed as (see Appendix B) 
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Figure 5.7 — Concentrated -forces on crack surface. 
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and 


T r0 ~ 
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The a rr and r r g from equations (23) and (24) were equated ,to 

A * 

corresponding sum of the stresses from equation (22) at "discrete points 
on the circular arc and coefficients A n and were determined from 
the resulting set of equations. Thus by superposition, the two sets of 
radial and tangential stresses cancelled each other at intervals along 
the circular arc of unit radius to effectively create a stress-free 
boundary. This superposition is illustrated in figure 5*9* The cr rr 
and r r Q on the unit circle in figure 5* 9(b), due to the concentrated 
loads, opposed the cr and T r0 in figure 5* 9(c) from equations (23) 
and (24), and, as a result, the stress-free hole with cracks in 
figure 3.9(a) was simulated. 

The stress-intensity factor for the case presented in figure 5* 6(c) 
and again in figure 5.9(a) was found by superposition of the stress- 
intensity factors for figures 5* 9(b) and 5* 9(c). The stress -intensity 
factor for each concentrated force was given by Paris Qjf] as 
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k 


2P \/a 

irt \/ a 2 - l 2 


(25) 


therefore the stress-intensity factor for figure 5.9(h) was 

M 

i=l 7Tb 

> 

where ( M' was the number of concentrated loads. The stress -intensity' 
factor- for the loading _ shown in figure 5. 9(c) was. written in terms of 
$(z) by -Paris and Sih (52J as 

k = 2\/2 lim V z “ a $(z) (27) 

z a 

Therefore from equations (26) and (27) we obtain,' 

M _ 

V 1 2P,* Ja _ • 

k = / . — - + 2 \/2 lim \J z - a $(z) (28) 

i=l irfc /a 2 - b? z a 

which represents the stress-intensity factor for the case shown in 
figures 5- 9(a) and 5* 6(c) and consequently for the original case, 
illustrated in figures 5«6(a) and 5*1* 

Equation (28) was evaluated with the aid of a digital computer for 
discrete crack lengths and for the corresponding values of cr(x). The 
stress-intensity factor was plotted for each crack length and a curve 
was drawn through these points to obtain a curve similar to that 
presented in figure 5.5 but for a general stress distribution, cr(x), 
along the crack plane. For an evaluation of this general approach 
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ct(x) was specialized for the elastic distribution and the stress- 

• / 

intensity results were compared with the curve from figure 5* 5* The 
close correlation in figure 5*10 between these results from equation (28) 
and Bowie's solution demonstrated the adequacy of the superposition and 
boundary collocation techniques applied in the development of 
equation (28). 

These results, also presented in table VI, show that beyond a 
crack length of 0.7 inch less than 6 percent difference exists between 
the stress -intensity factors for the two cases. Consistent with these 
results, for all loading conditions nearly identical crack-growth 
behavior was observed for cracks larger than 0.7 inch, as shown 
previously in figure 5.4. 

5.3.2 Calculation of Crack-Growth Curves From Stress -Intensity Solutions 

The correspondence' between stress -intensity factor and crack-growth 
rate demonstrated "by Figge and Newman [50] was utilized to calculate 
crack propagation curves. The basic crack-growth, behavior of 2024-T3 
aluminum alloy was experimentally established through laboratory tests 
of centrally cracked sheet specimens for which crack-growth rates were 
measured and the corresponding stress-intensity factors calculated 
from equation (20) using the Irwin-Westergaard finite width correction 
from (52J. Experimental results are presented in figure 5 . 11 to 
characterize the crack-growth behavior for this material. This curve 
in figure 5*11 together with the stress-intensity .solution for a 
given cr(x) were used to calculate crack length against cycles curves 
as follows,.- Stress -intensity factors were found for crack lengths a 
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TABLE VI. - CALCULATED- STRESS -INTENSITY FACTORS 










Crack props 



Figure 5.11 — Master crack propagation curve for 202k^S3 aluminum alloy, 



no 


and a + Aa and corresponding crack-growth rates were found from the 
master ’Curve in figure 5*H* The number of cycles required to extend 
the crack by Aa was determined from 


AN = 


2Aa 


da 

air 


da 

m 


a+Aa 


(29) 


The total number of cycles required to grow the crack to a given length 
was found by dividing the crack into intervals Aa, and by summing the 
AN ' from equation (29) associated with each Aa. For the present study, 
a digital computer was used in this procedure to calculate curves for 
crack length against cycles. These results are presented in the next 
section .and are compared with observed results of Section 5*2. 

5., 4 Comparison of Calculated and Observed Crack Propagation Behavior 
Comparisons are made in this section between calculated and 
observed crack-growth’ results for the reference case of constant- 
amplitude, S = 20 ksi, loading and for the case of two-level beneficial 

• i 

prior loading." The stress distributions, cr(x), on the transverse axis 
were found, respectively, by elastic and elastoplastic analyses for 
these two cases, and the methods from Section 5* 3*1 were used to 
develop, stress-intensity factor solutions. These stress-intensity 
results were then utilized together' with the master crack propagation 
curve in figure 5.11 to calculate crack length against cycles curves 
similar to the experimental curves in figure ^.b. 


in 


5*^*1 Constant-Amplitude (S = 20 ksi) Loading 

A system of concentrated forces was developed to simulate the 
elastic stress distribution on the transverse axis by the procedure 
presented -in Appendix C. This force system was used to develop the 
.stress -intensity curve previously presented in figure 5.10 for an 
evaluation of the accuracy of the stress-intensity analysis. This 
stress-intensity solution led to the crack length against cycles 
curve presented in figure 5.12. Because of the inherent scatter in 
crack propagation . data in figure 5* If upper and lower limits were 
used to characterize this crack propagation behavior as shown. These 
limits were used to calculate the corresponding range of crack-growth 
behavior for the case of constant -amplitude, S = 20 ksi loading. 

This expected range of behavior is indicated by the shaded zone in 
figure 5.12. The solid curve within this shaded zone represents 
average behavior found from the corresponding average curve in 
figure 5-11. 

The symbols in figure 5.12 represent test data. For the two tests 
conducted with no prior loading (constant-amplitude) one set of test 
results fall slightly outside of the range of calculated behavior. 

The second set of results is within this range 'and the average of 
these two sets of data was within the shaded zone. Because of large 
errors in crack length against cycles behavior that may result from 
rather small errors in stress -intensity factors (and therefore crack 
propagation rates) for short crack lengths, the correlation between 
experimental and calculated results in figure 5*12 is. exceptionally 


close. 1 



Normalized crack length, L Ip 



Figure 5;12 — Comparison of calculated and observed crack growth curves 
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5.4.2 Beneficial Two-Level Loading 

An approximate analysis is presented in' Appendix C for the stress 

5 ; ' 

distribution corresponding to beneficial two-level loading. This stress 
distribution was discretized and also simulated by a system of con- 
centrated forces. This force system together with the procedures 
outlined in Section 5* 3*1 resulted in the stress-intensity solution 
listed in table VI and plotted in figure 5*10« The low stresses cal- 
culated in Appendix C for small values of x produced stress -intensity 
factors that were lower than those for the constant- amplitude case. 

As expected, these low values of stress intensity created a large 
crack propagation life, as shown, in figure 5*12. Again the results are 
presented as a range of calculated behavior. For this case, however, 
the single experimental crack length against cycles curve did not lie 
within the calculated range. This discrepancy may have been the result 
of scatter in material behavior (note that data- points fall outside of 
the limits used in figure 5-U to characterize the assumed scatter). 
Additional tests are required for a critical assessment of this 
analytical approach, however, the single experimental curve in 
figure 5.12 illustrates reasonable- agreement between calculated and 
observed results. As previously mentioned, small variations in stress- 
intensity values for short cracks produce significant changes in 
predicted crack propagation lives. For this reason, minor refinements 
in the elastoplastic analysis of the stress distribution or in the 
simulation of this stress' distribution, by concentrated forces could 
noticeably improve the correlation between experimental and calculated 
results in figure 5*12., « . , 1 . ( 


6. CONCLUDING REMARKS 


An analysis of the role of a typical stress concentration in 
structural fatigue has -been presented. The analysis considers the 
cyclic stress-strain behavior at stress-concentration sites,, the 
initiation of fatigue cracks, and the initial stage of propagation 

for cracks growing from the stress-concentration site. The stress 

* 

concentration dealt with in this study was that produced by 1 a circular 
hole in a sheet specimen under remote load. The behavior of this' 
specimen, was investigated for constant-amplitude loading and two-level 
loading, that created stress interactions typical of those experienced- 
by structures in service. 

For the analysis of stress-strain conditions at the stress- 
concentration site, attention was focused on the point of maximum 
stress concentration, and conditions at this point were referred to 
as "local" conditions. Cyclic local stress and strain were .simulated 
in simple unnotched specimens by loading these specimens* between cal- 
culated limits for local behavior. This simulation, was performed on 

' < 

i , 

a cycle-by-cycle basis and,* as a result, residual stress effects and 

cyclic variations in material behavior., that generally complicate 

local conditions, were automatically accounted for in quantitative 

estimates of local’ stress arid strain; For -an evaluation of the 

simulation procedure, local cyclic stress-strain behavior was experi- 
■ . * ; * 
mentally determined 'by ( a; companion- specimen method and the two sets 

i • ' 

of results were compared. . 
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The simulation procedure was continued until the unnotched 
specimens failed by fatigue. Because the stress- strain conditions in 
these unnotched specimens were assumed to be identical to the con- 
ditions at the stress-concentration site, these fatigue lives were 
used as estimates for fatigue crack- initiation lives for notched 
specimens. Notched specimens were tested in- fatigue to obtain observed 
crack- initiation lives which, were 'used for quantitative assessment of 
stress-interaction effects predicted by the simulation procedure. 

The effects of prior high loading on cracks growing from stress 
concentrations were demonstrated by tests with notched specimens. ' 

Crack length against cycles curves were experimentally determined 
from these -tests and were compared with similar predicted curves. 
Predicted curves were calculated by expressing the state of stress 
near the crack tip in terms of stress-intensity factors and by using 
these stress-intensity factors .together with an experimental curve 
generally characterizing the crack propagation behavior of the specimen 
material. ' } . * < 


Based on analytical and experimental 'results presented in this 
study, the following conclusions are presented. 

*« t 

1. Prior high loading on a notched specimen may dramatically 
* \ . *\ < * ’ ‘ 

« * , * r 

alter the, local mean stress corresponding to subsequent loading. Prior 


loading may be .selected to 

r ' < 

stress. 


either ' increase or decrease this local mean 


i 

2. The companion -specimen method together with the strain- 
coupled servo-control system presented in this study provide means 
for experimentally investigating local cyclic stress and strain. 


3* The simulation procedure for approximating local cyclic 
stress- strain behavior adequately accounts for local plasticity. 
Furthermore, changes in stress-strain properties due to cycling are 
inherently incorporated in the calculations of local behavior. 

4. The increase in crack- initiation lives for tensile prior 
loading on notched specimens is the result of the net effect of 
beneficial compressive residual stresses and the detrimental effect 
of damage acceleration caused by large excursions of local strain at 
the stress-concentration site during prior loading. For compressive 
prior loading, detrimental tensile residual stresses and damage 
acceleration combine to decrease crack- initiation lives. 

3 . The simulation procedure for estimating crack-initiation 
lives for stress concentrations using unnotched control specimens 
accounted for the effects of stress interactions and, as a result, 
yielded reasonable estimates. These estimates were, however, conserva- 
tive for both constant-amplitude and two-level loading due to a size 
effect not accounted for in this simulation procedure. 

6. The initial stage' of crack propagation for cracks emanating 

from a stress-concentration site may be altered by prior loading during 

the crack- initiation phase. Prior tensile loading decreases initial 

crack propagation rates and compressive prior loading increases 

propagation rates. 

» 

7 . "Reasonable predictions of ; crack length against cycles were 
obtained' for the initial stage of crack growth when the stress state 

1 

near the crack tips were defined by stress-intensity factors which 
were used to establish crack propagation rates. 


117 


Although the procedures and conclusions in this study are applicable 
only for sheet specimens -with stress concentrations, it is possible that 
they will prove useful in fatigue analysis of more complicated con- 
figurations and loading. 
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8. APPENDIX A - DEVELOPMENT OF STRESS EQUATIONS 
FOR CONCENTRATED FORCES ON A CRACK BOUNDARY 


Irwin [55] presented the Westergaard stress function 

_l/2 

2P^ fo ,0 

Z l(z) = 


:rtt(z 2 - b 2 ) [z 2 - a 2 J 


‘a 2 - h 2 


(Al) 


for the problem shorn in figure 5.7. The- stress components in 
Cartesian coordinates were written as 


' a xx = ReZ l " 

Cyy' = ReZ x + ylmZp 


; 


and 


V = -y ReZ i 


j 


For the purpose of calculating stresses along the circle of unit 
radius equation (A2)' was transformed to polar coordinates by 


cr rr = cos 2 0a xx + 2 sin 9 cos + sin 2 0a^ 


(A2) 


'N 


and 


) (A3) 


T r9 = -sin 0 cos + (cos 2 0 - sin 2 9)T xy . + sin 0 cos 00^ 


The substitution of -equation (A2:) into equation ,(A3) yielded 
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(ReZ^ - yl^ZjJcos 2 © - 2yReZj[ sin 

+ \kez 1 + yi ril zi]sin 2 e 


y[ 2 Vi sin 8 cos 9 - ReZ-j,(cos 2 0 


8 cos 0 


- sin 2 9)J 


9. APPENDIX, B - DEVELOPMENT OF STRESS EQUATIONS 
i'OR BOUNDARY COLLOCATION PROCEDURE • 


Muskhelishvili |j>5j ' expressed, stress components in Cartesian 
coordinates in terms of the stress functions 0(z) ‘and f(z) by 
the equations . t 

\ 4 r 


a ll + a 22 = 2 &' ^ + ( z ) 3 


and 


’ > . . ‘(EL) 




^22 ” CT 11 + 2 i T i 2 *=; 2 '^"(z : ) + r ('z)] 
or with ' • • t • 1 ‘ . 

■i ‘ j 

® (z) = (z), $(z)-= (z) * 

I * 

Muskhelishvili rewrote equations (Bl) as •> 

a ll + a 22 = 2 .&( z ) + $ ( z )3 

and 

a 22 ” a U + 2iT i2 = (z) + $(zj] 

The introduction of D(z.) also from 


* i 






j 


D(,z) = $(z) + z$ f (z) + $(z) 


(B2) 


(B3) 


(B4) 


yields 
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a ll + a 22 = 2 GK' z ) + ®(z)] 




and 




(B5)' 


a 22 ~ ff ll + 2iT 12 = 2 B Z - z) $t (z) - ®( z ) + 0 (b)] 

In cylindrical coordinates equations (B5.) become 
*rr + °08 = 2 &( z ) + 

and ^ (B6) 

age - a rr + 2iT rQ = 2 Q(z - z)$* (z) - $(z)- + fl(z)]e 210 
Subtracting equations (b 6)- yields 


cr rr - - i T r e = ®(z). + *(z) - [(z - z)«* (z). - #(z)_+ fi(z)]e 2ie (B7) 

Newman presented expressions for $(z)*,' and- fl(z) in reference, 
following a general form developed by Erdogan for a plate con- 
taining a, crack. From QjiT] \ 




•M 

®(z) = ^ 


■n 


V< 


n= p (z 2n "^\/z 2 - a 2 z 211 ^ ' 


and 


(b8) 


MV' 1 " 

V^I I . * < 

0(b) = Y {— r- = r 

n=l\ z2n '"V z2 ~ a 


^2n 
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where a is the half crack length and An and Eft are arbitrary 
constants. For the present study these Muskhelishvili stress functions 
were substituted into equations (B7) and real and imaginary parts were 
separated to obtain expressions for cr rr and T r g. 

Equations (b 8,) were rewritten in exponential form, using a change 
of variables. 


z 

■z 

z 


+ 


re-*-® 
a - r^e 
a = rge 





‘ (B9) 


as illustrated in figure Bl. From equations (B8) and equations (B9) 
one -obtains 


M 


$(z.) - ^ 


A n exp \ -i 


.(2n - 1)0 +4p_+- ^ 


An e 


-2in9 


n=l 


T . r 


2n-l 


\/W2 


_2n 


(ELO) 


Differentiating $(z) from equations (b 8)' and substituting from 
equations (B9) yields ' , 


M 


(z) = 


n=l 


A n (-2n + !)■ 


exp 


r2n \/ x l r 2 ‘ 




A n exp \ -i 


(2n - 2)8 + + ^| 2 


2n-2 




■ “• e ^' (- i ( 2n + 1 ) 9 > 


(ELI) 
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The conjugate <±>(z ) is 


M 


$<z) = 


n=lj 


Ay 


„2n-l 


2 


exp 


(2n - 1)6 * ~2 * ~2 


M 

r 


2±n8 

2n 


(B12.) 


and 


£2(z) = fi(z) 


(HL5) 


By substitution from equations (b 8), (BIO}, (ELL), (B12), and (B13) 
into equation (B7) with the introduction of p-^ and p 2 






and 


) • ,(B14) 




J 


one obtains 




V 
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M 


- ir„ Q = ^ A 
n=l 


rr re / . 




exp < -ijj2n -1)0+0 


lJ'l 


+ exp /i [(2n -1)0 + p j\ . - 2l '^ g n ~ l) exp. f-i [(2n - 2) 6 *+0 

V. J r^ny r ^ r ^- - v. 


2iy 


fin-2 


M 


\ Z ( r ! r 2-) 5 


exp [~(2n. --4),0 +p 2 


I *■{*[•-“ 

n=l 


-2in0 + e 2in9^ , ^.ny eX p f -i;(-2n -1)0^ 
r 2n+l v 


+ “In exp {" i '( 2n - 2 ) e } 


-(•EL5) 


. * 


Separation of' real and imaginary parts of equation *(B15) leads to the 
following trigometric forms- for - a xr and T r9* 


M r f . • 

'""In 8 ** ^ °° s£(g11 '-' 1 




i 1 

v * 


r “/ r r r 2 


l ‘ 


y 


2n-2’ 


\/Crir 2 )3 


sin [(2n -4)0 + 0 p~) ; 


* 2 ®n 


'£ 
• r 2n 


cos 2n0 


y sin (2n-l)0 + cos(2n-2)0 


2n+l 


2n 


(Bl6) 


and 
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T r0 “ 


m r 

£{”* 


^ COE gsn - 2)0 + pj . 


r 2 "/ri.r 2 


y 


„2n-2 


>/W 


cos [j2n - 4)0 + pj 


+ 2B, 


n 


^2- cos(2n - 


„2n+l 


+ _i~ sin(2n - 2)6 
r 2n ' 


Equations (Bl 6 ) and (EL 7 ) axe the expressions for a rr and T r g 
used in the "boundary collocation procedure. 



1)6 

(B17) 



10. APPENDIX e - CALCULATION OF CONCENTRATED FORCE SYSTEMS 
SIMULATING STRESS DISTRIBUTIONS ON TRANSVERSE 
AXIS OF SPECIMEN 

10.1 Elastic Distribution for S = 20 ksi 

The elastic stress distribution for the specimen configuration 
used in this study was nearly identical to that for an infinite sheet 
•as. illustrated pre-viously in figure 3*2. Therefore, for convenience 
the elastic solution Q>8] 



for an infinite sheet with a circular hole of radius p was used. • 
Equation ' (Cl) was. evaluated for a = 20. ksi and . 0 ■=• 0, and Oqq 
was plotted for the full width of the specimen.* _ As typically 
illustrated in figure '.Cl, this stress distribution was divided into 

- t , * * 

discrete segments. The avei-age, stress (ct 0 0 for each segment was 

• • ' ' t 1 , * 

used to calculate -a statically; equivalent concentrated force . from 

P ± = Abt^aee).^ ^ ■ ( C2 ) 

where t is specimen thickness and £b = 0..0025 inch. By this 
procedure a system of concentrated forces was developed to .approximate 
the elastic stress distribution for S = ,20 ksi. 

10.2' Elastoplastic Distribution For Beneficial Prior Loading 

The prior loading produced localized plasticity as previously 
illustrated in figure k.k. The procedure used to calculate the stress 
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distribution in figure 4.4 is described in this section. The elastic 
stresses corresponding to the prior loading stress level (S = 40 ksi) 
were determined from equation (C1-) and 


J rr r 


Ug 

2 


/ \ 2 

/_ \2 / \4 


& ) - 

1 - 4(£) + 3[£) 

cos 20' 

W 

W W 


I 

.. _ 



(C3) 


along the transverse axis. The effective stress o a corresponding 

e G '. 


to this elastic solution was found from- 


+ % 


CT eo = [ a rr D “ a rr o a 0e o + a ee] 


l/2 




•An approximate elastoplastic analysis was based on’ these elastic stresses 

1 

using the approach presented by. Dixon and Strannigan^ { 25 ] • • From 

J 1 1 

equation (ll) in Chapter 2 an expression for. ctqq was written as 


°ee - 0 ee o y E 



(C5) 


where E s is the secant modulus from the uniaxial stress- strain curve 
corresponding to the effective stress a e . This effective stress o e 
was found from 


cr P a, 




eo 



E 


(c 6.) 


The elastic effective stress a eo was found from equation (C4), but 
because E g depended on’ a e , equation (C 6 ) could not be solved 
directly for the E s required in equation (C 5 ). A graphical procedure 


J 


135 


was used to obtain E s , as illustrated in figure C2. The curve of a 
against yE s /E was constructed from the uniaxial stress-strain curve 
for 2024-T3, and equation (c6) was plotted for a given cr eo . The 
intersection of these two curves represents a simultaneous graphical 
solution for equation (C6) and the stress- strain relationship. This 
value of \J E s /E and ctqq. , corresponding to the a eo used in 
figure 02, were substituted into equation (C5) to obtain <j 00 . To 
account for the cyclic strain hardening displayed by 2024-T3 ‘(see 

t 

for example fig. 3* 20 ), the cyclic stress- strain curve from figure 3*16 

was used in this graphical procedure. f 

* 

The above procedure was repeatedly applied to- -construct the 

. • k 

plastic portion of the stress distribution in figure C3 for 
S = 40 ksi. Because plasticity ''was accounted for by .attenuating 

' ' .. i ’* ‘ 

elastic stresses this procedure made no correction to the elastic 
portion of the elastoplastie stress distribution. To satisfy equili- 

1 ’ < t < * 

brium the elastic portion of- this stress distribution must be altered 
to account for the decrease in area under the curve caused by the 
attenuation of stresses due to plasticity.- For the present analysis 
the elastic portion of the curve was raised uniformly until the 
increase in area under this part of the curve was approximately equal 
to the above mentioned decrease' in area. This minor alteration to the 
stress distribution is shown in figure C3> 

The residual stress distribution in figure C3 "was found by 
subtracting stress range values corresponding to elastic unloading 
from the elastoplastie stress distribution, and by applying the above 





Figure C3 — Stress distribution on transverse axis for beneficial prior loading 


procedure to make a small correction for reversed yielding • Hie 
resultant .stress distribution (for application of S = 20 ksi, 

following the unloading from S =.'l0 Ksi) was obtained by adding 

1 ' ■, 

elastic stress range values tor- S.= 20 ksi to the residual stress 

if 

distribution. 

. , ^ 

. *■ ' i f 

i * < ’K 

The system of concentrated ' forces simulating. this stress 

< > * ’ f t 

' . * j 

distribution was calculated using equation (C2) and’.* the procedure' 

' ’ , - _i 


in figure Cl. 
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